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The  properties  of  oxide  layers  on  smooth  and  rough  aluminum 
surfaces  are  considered  in  this  study.  Using  ellipsometric  techniques, 
the  layer  thickness  and  the  four  parameters  describing  the  complex  in- 
l dices  of  refraction  for  film  and  substrate  are  assigned  effective 

values  and  interrelations  are  deduced.  The  organizing  feature  of  the 
work  is  the  use  of  data  trajectories  generated  by  determining  the 
I ellipsometric  parameters  ip  and  A over  and  over  for  samples  that  were 

heated  between  measurements.  The  heating  caused  increases  in  the  oxide 
layer  thickness,  so  that  the  data  trajectories  could  be  used  along  side 
I theoretical  thickness  curves  to  determine  appropriate  effective  media 

parameters.  One  of  the  more  interesting  results  was  the  representation 
of  rolling  grain  marks  on  alclad  by  an  effective  complex  index  of 


t refraction.  The  effective  index  determined  was  dependent  on  the  sample 

orientation. 
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SECTION  I 
INTRODUCTION 

When  aluminum  is  exposed  to  air,,  oxidation  occurs  almost 
instantly.  Thin  films  of  aluminum  oxide  form  on  all  surfaces  of  alu- 
minum. When  heated,  the  thickness  of  the  oxide  films  increases. 

Ellipsometry  is  a well  established  method  of  measuring  properties 
of  thin  films  on  optically  flat  surfaces.  The  use  of  ellipsometry  for 
rougher  surfaces  is  unusual  and  there  is  no  established  method  for 
determining  the  relevant  optical  properties  of  the  film  and  substrate. 
These  properties  must  be  known  in  order  to  determine  the  film  thickness. 

A new  method  for  determining  these  properties  is  the  use  of  a data 
trajectory  found  by  repeatedly  heating  and  then  measuring  the  values  of 
ip  and  A for  a given  sample.  The  thickness  (and  roughness)  of  an  oxide 
film  will  grow  with  heat  and  a series  of  data  points  is  generated, 
characterized  by  increasing  (unknown)  thickness.  This  method  was  first 
used  by  Cap.  M.  Lydon[l]  at  Texas  Tech  University  in  1975.  This  trajec- 
tory can  be  used  to  relate  the  growth  of  the  oxide  film  to  the  complex 
refractive  indices  of  film  and  substrate: 

nf  = nf(l  + jKf)  and  n$  = ns(l  + jK$) 

This  relationship  is  very  complicated  and  requires  the  use  of  a 
computer  for  the  calculations. 

The  study  of  optically  flat  aluminum  surfaces  was  undertaken  in 
this  work,  using  data  trajectories  as  described  above.  In  order  to 
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interpret  the  data,  a policy  of  using  the  data  trajectories  to  define 
functional  relationships  between  media  parameters  has  been  developed. 

It  was  found  that  the  oxide  roughness  on  nearly  ideal  aluminum  sur- 
faces can  be  modelled  by  taking  the  film  extinction  coefficient  K$  to 
be  a function  of  film  thickness,  t: 

and  by  imposing  a relation  between  the  real  and  imaginary  parts  of  the 
aluminum  substrate  index: 

Ks  = 7.175  - | ns 

This  relation  also  interpolates  all  widely  scattered  values  for  the 
index  of  aluminum  found  in  the  literature.  The  actual  thickness  of  the 
oxide  film  is  the  key  to  picking  the  correct  values  of  n$  and  K$.  As 
the  data  collected  in  this  study  was  compared  with  work  by  others,  a 
reasonable  choice  for  the  aluminum  refractive  index  was  selected: 

n$  = 1.08(1  + j 5 . 555 ) 

Application  of  the  trajectory  idea  to  the  rougher  surface  of 
alclad  aluminum  indicated  that  rough  surfaces  could  also  be  modelled  by 
variation  of  index  parameters.  It  was  found  that  substrate  surface 
roughness  could  be  modelled  by  increasing  the  value  of  ns  and  requiring 
to  be  a function  of  n$: 

Ks  = 9.75  - 4.77ns  + 0.505n| 
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Also,  the  effective  value  for  was  found  to  be  a function  of  both  t 
and  ns  when  the  substrate  surface  is  rough: 

Kf  ■ (°-°'>Wns 

The  ideas  employed  in  the  present  work,  using  data  trajectories 
to  determine  effective  parameter  relationships,  show  great  promise  for 
work  with  thin  films  on  both  smooth  and  slightly  rough  surfaces.  For 
example,  it  was  possible  to  associate  the  present  roughness  criterion 
to  the  orientation  of  rolling  grain  marks  on  the  alclad  samples.  This 
interesting  check  is  described  in  Section  VI. 

In  Section  II  the  reflection  coefficients  and  their  relation  to 
the  film  and  substrate  parameters  n^,  Kf,  n$,  K$  and  t are  derived. 
These  are  then  related  to  the  parameters  ip  and  A.  The  mechanics  of 
the  actual  ellipsometric  measurement  of  the  values  of  ip  and  A are 
discussed  in  Section  III.  In  Section  IV  an  analysis  of  errors  in  the 
ellipsometric  measurements  and  a sensitivity  study  for  the  various 
parameters  are  given.  The  study  of  optically  flat  aluminum  surfaces 
is  described  in  Section  V.  Finally,  the  study  of  the  oxide  layers  on 
alclad  alloy  is  described  in  Section  VI. 


SECTION  II 


THEORETICAL  BASIS  OF  ELLIPSOMETRY 

Ellipsometry  is  a method  of  measuring  the  changes  in  light  as  it 
is  reflected  by  a surface  in  order  to  determine  the  physical  properties 
of  that  surface.  To  understand  what  changes  as  light  reflects  from  a 
surface  and  how  this  is  related  to  that  surface,  it  is  necessary  to 
study  some  of  the  basic  properties  of  light.  This  chapter  gives  a 
brief  description  of  the  laws  governing  the  reflection  and  refraction 
of  light. 

In  ellipsometry  the  light  used  is  a plane  wave  of  monochromatic 
light.  The  surface  on  which  measurements  are  made  must  be  smooth  and 
flat  with  respect  to  the  wavelength  of  light.  It  must  be  remembered 
that  measurements  are  made  with  light  at  only  one  wavelength  and  the 
index  of  refraction  can  be  different  at  other  wavelengths. 

Light  is  part  of  the  electromagnetic  spectrum  and  is  governed  by 
Maxwell's  Equations  which  describe  the  electromagnetic  field.  These 
equations  relate  the  electric  vector,  ?,the  magnetic  induction,  ft,  the 
electric  displacement,  ft,  the  magnetic  vector,  ft,  and  the  electric  cur- 
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These  five  basic  quantities  are  further  related  by 


5 = o t 

(2.5) 

ff  = et 

(2.6) 

t = ptf 

(2.7) 

Here  a is  called  the  conductivity,  e is  known  as  the  permittivity,  and 
p is  called  the  permeability.  These  three  terms  describe  the  proper- 
ties of  the  media  in  which  the  electric  field  is  propagating.  When 
o=0  the  media  is  non-conducting  and  the  electric  conduction  current 
density,  j,  vanishes.  The  electromagnetic  components  are  not  damped  as 
the  wave  propagates  through  the  media  and  the  material  appears  trans- 
parent. The  properties  of  light  are  easier  to  describe  for  nonabsorb- 
ing (non-conducting)  media  than  for  a conducting  material.  The  laws  of 
reflection  and  refraction  will  first  be  described  for  surfaces  with 
o = 0. 

For  a medium  that  is  homogeneous  and  contains  no  charges  or 
currents,  j = 0 and  o = 0,  Maxwell's  Equations  can  be  used  to  derive 
two  wave  equations: 

V2  t - ep  |^|  = 0 (2.8) 

and 

V2  ft  - ep  = 0 . (2.9) 

When  the  field  is  monochromatic  wi  th  angular  frequency  w,  the  time 
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derivative  = - iw  for  ? and  It  of  the  form  £ = E0e~iajt.  The 
solutions  for  the  above  wave  equations  are 

t - E0ei(k2'  - (2.10) 

and 

I * = H0el(kz'  ' ut)  (2.11) 

where  the  propagation  is  in  the  z'  direction  and  k = >/c\i  w. 

When  light  falls  upon  a boundary  between  two  media  and  is 
transmitted  into  the  second  media,  the  wavefront  is  continuous,  al- 
though it  has  a kink  at  the  boundary  as  shown  in  Fig.  (2.1).  The  wave 
appears  to  be  bent  at  the  boundary.  The  velocity  of  propagation  of  the 
wave  changes  at  the  boundary  and  the  amount  of  refraction  or  bend  is 
related  to  the  change  in  velocity.  The  ratio  of  the  first  velocity  to 
the  second  is  called  the  refractive  index  for  refraction  from  the  first 
media  to  the  second. 

ni2  = ^ (2.12) 

v2 

An  "absolute  refractive  index"  can  also  be  found,  it  is  the  ratio  of 
the  velocity  of  light  in  a vacuum  to  the  velocity  in  the  media, 

n * ~ . (2.13) 

If  the  absolute  refractive  indices  of  two  media,  nj  and  n2  are  known, 
then  the  refraction  from  the  first  media  to  the  second  can  be  found. 


(2.14) 


The  refractive  index  or  the  index  of  refraction,  as  it  is  often  called. 


is  defined  to  be 


n=  jSE 

VrAUn 


(2.15) 


where  e0  and  y0  are  vacuum  permittivity  and  permeability. 

The  laws  of  reflection  and  refraction  can  be  derived  from  Maxwell's 
Equations,  but  they  were  discovered  experimentally  long  before  Maxwell's 
time.  The  incident  beam  is  assumed  to  be  a plane  wave,  the  wavefronts 
being  normal  to  the  incident  ray.  The  reflected  and  refracted  rays  lie 
in  the  plane  formed  by  the  incident  ray  and  the  normal  to  the  surface  at 
the  point  of  incidence.  This  plane  is  shown  in  Fig.  (2.2).  The  angle 
of  incidence  (6^),  the  angle  of  reflection  (0r),  and  the  angle  of  re- 
fraction (6t)  are  measured  between  the  normal  to  the  surface  and  the 
appropriate  ray.  For  reflection: 


9i  = 9r 


(2.16) 


For  refraction: 


m sine^  = n2  sin©t 


(2.17) 


This  is  often  called  Snell's  law. 

To  fully  describe  reflected  and  refracted  waves  it  is  necessary  to 
know  more  than  the  angle  of  propagation.  The  amplitude  of  these  waves 
must  also  be  found.  Expressions  describing  the  amplitudes  of  the  re- 
flected and  refracted  waves  are  found  by  imposinq  the  boundary  condi- 
tions. These  conditions  demand  that  the  tangential  components  of  £ and 
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H must  be  continuous  across  the  boundary. 

For  convenience  in  discussing  the  boundary  conditions,  the  electric 
and  magnetic  fields  will  be  described  in  terms  of  amplitudes  ?{t)  and 
ft(t),  obtained  by  extracting  the  propagation  factor  exp(ikz'): 


t = ?(t)e1kz' 
B » /)(t)e1k2' 


(2.18) 


where  z'  measures  distance  in  the  direction  of  propagation.  As  shown  in 

Fig.  (2.2),  the  incident  plane  wave  propagates  in  the  direction  of  the 

~(i)  ~(r)  ~(t) 

unit  vector  ev  . The  unit  vectors  ev  and  ev  ' represent  the  direc- 
tion of  propagation  of  the  reflected  wave  and  the  transmitted  wave, 
respectively. 

The  tangential  field  components  can  be  described  in  terms  of  the 
coordinate  system  shown  in  Fig.  (2.2).  The  z axis  is  normal  to  the 
boundary,  the  x axis  is  on  the  boundary  and  in  the  plane  of  incidence, 
and  the  y axis  is  on  the  boundary  but  normal  to  the  plane  of  incidence. 
In  order  for  the  tangential  components  of  t and  ft  to  be  continuous,  the 
following  conditions  must  be  satisfied: 


H (i) 

x 

H <*> 


+ E (r)  = 


t Ey(r)  . 

+ H <r) 


= H 


♦ « (r) 

y 


= H 


(t) 

(t) 

(t) 

(t) 


(2.19) 


> 


I 


> 


* 


» 


I 


» 


I 
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For  the  purposes  of  ellipsometry,  the  components  shown  above  are 
conveniently  replaced  by  related  components,  which  will  now  be  intro- 
duced. The  y-components,  perpendicular  to  the  plane  of  incidence,  will 
be  called  the  s-components,  a trivial  change  of  name.  The  remaining 
components  of  the  amplitudes  of  the  three  beams,  however,  will  be 
described  in  terms  of  three  unit  vectors  p^  where  a represents  i,  r 
or  t: 


p(r) 

p(t) 


= y x ^ = s X 
= J X elr)  = s X c(r) 


= 9 X c(t>  = s X e(t) 


(2.20) 


where  s = y. 

Using  these  unit  vectors,  the  electric  field  amplitudes  of  the 
three  beams  can  be  written: 


+ A 

£(r>  = R$s  + Rpp(r>  (2.21) 

t(t)  • Tss  ♦ yW 


Now,  in  the  x,  y,  z coordinate  system  the  components  of  the  three 
directions  of  propagation  e and  the  three  related  vectors  Cp^  are: 


/ sine.  \ 


:(i)  s 


\ 


sine 

I 

0 

-cose 


;{i>= 


/-  cose^ 

o 

.Sine., 

\ / x,y,z 


x,y,z 


i 


:(r)  = 


sine^ 


cosei 


x.y.z 


cosei 


. - o 


-sine.- 


x.y.z 


(2.22) 


~(t)  . 


x.y.z 


x.y.z 


Substituting  Eqs.  (2.22)  into  Eqs.(2.21),  one  obtains: 


*«>  ■ 


-ApCOSO, 


-A_sin0. 
^ P i 


x.y.z 


|(r)  = 


Rpcosei 


-Rpsin9ij 


x.y.z 


(2.23) 


*(t)  = 


r_Tpcoset 


-T^sine. 

P t. 


x.y.z 


The  magnetic  fields  and  may  be  obtained  from  the 


relation: 


(2.24) 


) 
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where  a represents  i,  r,  or  t.  Using  Eqs.  (2.22)  and  (2.23)  in 
Eq.  (2.24)  one  obtains: 


a«> . /s 

v Pi 


x,y,z 


tf(r)  = /n 

V Pi 


nW  = /S 

V U2 


fV'A 


/-Rscosel\ 


\ Rssin0i  / 
v / x,y, 

/ Tscos0t\ 

TP 

\Tssinet  / 

\ / x,y. 


/h  ^ 

x 

V" 

. (D 


I Hx<r)\ 

H (r) 

y 

V z / 

/Hx(t) 

M (t) 


(2.25) 


V 


(t) 


The  representations  in  Eqs.  (2.23)  and  (2.25)  may  now  be  substi- 
tuted into  boundary  conditions,  Eq.  (2.19)  to  obtain: 


[-Ap  + Rp]cosei  = -Tpcos9t 


As  + Rs  ■ Ts 


n^As  - R$]cosei  = n2Tscoset 


(2.26) 


n i [Ap  + Rp]  = nzTp 


where  the  approximations 


Pi  » P2  {5  Po 


r 
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have  been  made,  so  that 


/n  / /n  = /nn  = /£  1^1  = dj 

V yi  / V U2  V c2uq  V e2y2  n: 


(2.27) 


From  these  Eqs.  (2-27),  it  can  be  observed  that  the  parallel  and 
perpendicular  components  of  the  lightwaves  are  independent  of  each 
other.  These  four  equations  can  be  solved  for  the  amplitude  of  the 
reflected  and  transmitted  waves  independently  as  a function  of  the 
incident  wave  amplitude  for  both  the  parallel  and  perpendicular  cases: 

T 2 in!  cos  0. 

t = -H  = 1 

p " Ap  n2  cos  6i  + ni  cos  6t 


. _ s 2 ni  cos  0. 

Ls  = a = ! 

s ni  cos  0..  + n2  cos  0t 


Rp  n2  cos  0^  - ni  cos  ©t 

rp  “ Ap  n2  cos  0.j  + ni  cos  0^ 


Rs  ni  cos  0.  - n2  cos  0t 
rs  " As  ~ ni  cos  0^  + n2  cos  ©t 


(2.28) 


The  above  equations  are  called  the  Fresnel  Equations  and  the  quantities 
rp,  rs>  tp,  and  ts  are  called  Fresnel  reflection  coefficients  and 
Fresnel  transmission  coefficients. 

When  the  conductivity  a is  not  zero  these  equations  are  changed 


slightly.  The  wave  equation  is  changed  to 

n2  ? - ■ 9? 


(2.29) 
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The  plane  wave  solution  is 


t = t e1(k  z'  - ut) 


(2.30) 


where 


k2  = o)2y  (e  + i -j-) 

Cl) 


(2.31) 


All  of  the  previous  equations  are  the  same  for  a conducting 
medium  if  the  dielectric  constant  e is  replaced  by  a complex  term 


e = e + i — 
w 


(2.32) 


In  addition,  the  phase  velocity,  v,  and  the  index  of  refraction,  n, 
also  become  complex: 


v s£.  = _L 


(2.33) 


rUE 


and 


n = 


= /jE  = + iiC 

V'yoE0  / li o £ o 


VoEo  ,v  PoEo 

The  complex  index  of  refraction  can  be  divided  into  a real  and 
imaginary  part, 

n = n(l  + iK) 


(2.34) 


(2.35) 


here  K refered  to  as  the  extinction  coefficient  [2],  [3],  The  two 
real  quantities,  n and  K,  defined  in  Eq.  (2.35)  can  be  related  to  the 
media  parameters  by  equating  expressions  for  n2  obtained  from 


» 


Eq.  (2.34)  and  Eq.  (2.35): 


(1 

yo^o 


= n2(l  - K2  + 2iK) 


Comparing  real  and  imaginary  parts,  one  finds 


n!(l  - K2)  = rrf- 

PO^O 


2i/  , oy  _ c oy 
2e0yow  2w 


Solving  Eqs.  (2.36)  for  n and  K,  one  finds: 


* * *r  wtsf  - 13  - = 


n = c / oy  ~ c /us 
Y 2uK 


(2.36) 


(2.37) 


where  the  approximate  forms  are  useful  when  K«1 ; i.e.,  when  (o/eu)«l. 

Using  the  notation  introduced  above,  the  complex  propagation  con- 
stant ir  Eq.  (2.30)  can  be  expressed  in  the  form: 


O)  _ _U)~  _ 


(1  + 1 K) 


(2.38) 


Thus,  the  plane  wave  solution  of  Eq.  (2.30)  can  be  written: 

- ^Kz  i(^nz  - tot) 
t - f.  e c e c 


(2.39) 
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Tz 

The  first  exponential  term  in  this  expression  is  of  the  form  e where 
the  power  damping  constant  is  given  by 

a = = y^nK  , (2.40) 

c Ao 

with  A0  being  the  vacuum  wavelength  of  the  light.  For  metals  with 
large  conductivity,  K is  large,  so  that  the  field  is  damped  so  strongly 
that  it  can  penetrate  only  a fraction  of  a wavelength  into  the 
material . 

It  should  be  mentioned  that  the  expression  for  K given  in 
Eq.  (2.37)  has  a severe  defect  in  that  the  derivation  has  presumed  that 
(a/co))  ^ 0.  If  one  defines  this  ration  to  be  A: 


then,  according  to  Eq.  (2.37), 

Y 

so  that 


0 < K ^ 1 for  all  0 . (2.41 ) 

Nevertheless,  many  metals  which  are  good  conductors  must  be  described 
by  values  of  K that  are  considerably  larger  than  one,  say  K 
approximately  three  or  five.  This  point  is  discussed  at  some  length 
in  Born  and  Wolf [4].  It  is  actually  necessary  to  consider  e to 
be  negative.  This  unpleasant  feature  arises  because  the  media 
parameters  actually  depend  on  frequency.  In  other  words.  Maxwell's 
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Equations  are  not  truly  linear  in  the  fields  so  that  the  Fourier 
decomposition  shown  in  Eq.  (2.30)  is  not  strictly  appropriate.  The 
results  shown  in  Sections  V and  VI  of  this  study,  in  fact,  demonstrate 
values  of  K considerably  larger  than  one. 

The  effects  of  a complex  index  of  refraction  on  the  reflection 
and  refraction  of  light  at  a boundary,  are  such  that  the  basic  formulas 
all  have  been  exactly  the  same  form.  One  need  only  interpret  the 
various  parameters  as  being  complex.  The  law  of  refraction,  or  Snell's 
Law  becomes 


n2  sin  6t  = nj  sin 


(2.42) 


for  the  case  in  which  light  travels  from  a media  such  as  air  (o  = 0)  to 
some  media  where  the  index  of  refraction  is  complex  (o  f 0).  Since  n2 
is  complex,  the  angle,  0t,  must  also  be  complex  and  will  henceforth  be 
denoted  as  ©t.  The  Fresnel  Equations  are  also  effected  by  a complex 
index  of  refraction  as  well  as  by  the  complex  angle.  Looking  first  at 
the  reflection  coefficients, 

n2  cos  9i  - nj  cos  0t 


rp  “ n2  cos  ei  + nt  cos  St 


(2.43) 


ni  cos  0^  - n2  cos  0t 

r =• J L 

s ri!  cos  0.  + n2  cos  S. 


(2.44) 


Since  ri2  and  are  now  complex,  the  reflection  coefficients  and  ?s 
are  also  complex.  The  transmission  coefficients  also  become  complex. 


2 ni  cos  0i 


Lp  n2  cos  0i  + ni  cos  St 


(2.45) 
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2 n i cos  6. 

* = ' 

s ~ nj  cos  e,.  + n2  cos 


(2.46) 


The  complex  values  of  the  reflection  and  transmission  coefficients  mean 
that  the  ratio  of  transmitted  or  reflected  light  amplitude  to  the  inci- 
dent light  amplitude  is  complex.  This  means  that  a phase  change  occurs 
when  light  is  transmitted  or  reflected  from  a boundary  where  one  or 
both  of  the  media  is  complex.  Now  that  it  is  clear  that  the  effects  of 
conductivity  (o  t 0)  can  be  accoi-nted  for  by  simply  understanding  the 
indices  and  coefficients  to  be  complex,  the  notation  will  be  simplified 
by  dropping  the  wavy  bar  over  complex  quantities;  e.g.,  rp  e r . 

A full  description  of  all  properties  of  the  light  reflected  by  and 
transmitted  through  a media  interface  would  require  determination  of 
all  four  of  the  complex  Fresnel  coefficients.  This,  however,  is  more 
information  than  can  be  obtained  from  ellipsometry  alone.  The  ellip- 
sometry  procedure  described  in  the  next  section  is  able  to  determine 
only  the  complex  number  which  is  the  ratio  of  the  complex  reflection 
coefficients  rp  and  r$; 


r R A . 

= (tan^)elA 

rs  srtp 


(2.47) 


With  an  el lipsometer,  the  angles  \p  and  A can  be  determined  rather 
easily. 

The  behavior  of  light  when  it  reflects  from  a surface  with  a thin 
film  can  also  be  predicted.  The  boundaries  at  both  surfaces,  air-film 
and  film-substrate,  are  involved  in  this  case.  The  change  in  light  due 
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to  reflection  can  again  be  described  in  terms  of  angles  \p  and  A.  It 
should  be  noted  that  it  is  the  total  reflected  light  which  is  compared 
to  the  incident  light.  The  quantities  4/  and  A can  again  be  described 
in  termr  reflection  coefficients,  but  the  coefficients  come  from  two 
boundaries  and  the  terms  are  more  complicated. 

For  this  case  some  light  is  reflected  off  the  film,  but  some  light 
is  transmitted  through  the  film  and  reflected  from  the  third  media  and 
transmitted  back  through  the  film  as  illustrated  in  Fig.  (2.3).  The 
reflected  light  may  be  considered  to  be  made  up  of  a series  of  separate 
beams.  The  first,  reflected  at  the  first  boundary,  is  designated  as 
Rj  and  is  simply  found: 

Rx  = A r01  . (2.48) 

The  rest  of  the  beams  are  transmitted  through  the  film  and  reflected  at 
the  next  interface.  The  second  beam  is  the  first  one  transmitted  at 
the  thin  film  - air  interface  and  is  designated  R2.  The  other  beams 
are  reflected  between  the  two  interfaces  until  they  are  transmitted 
through  the  thin  film  - air  interface.  These  successive  beams, 

designated  as  R3,  R4 , R5, have  a smaller  and  smaller  intensity 

depending  on  the  reflection  coefficients: 

R2  = A t0,r12t10  e‘l26  (2.49) 

R,  - R toiri2riotio  e'i4!  • (2-50) 


» 


The  light  beams  transmitted  through  the  thin  film  are  delayed  in  phase 


Fig.  2.3  Reflection  of  Light  from  Two  Boundaries 
(Substrate  with  Film) 
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by  an  amount  26  for  each  pass  back  and  forth  through  the  thin  film 
because  of  the  extra  time  needed  to  travel  the  extra  distance.  This 
phase  delay  is  easily  determined  to  be: 


26  = 


/n(  - sin'1  0. 


(2.51) 


where  d is  the  film  thickness.  The  reflection  and  transmission  coef- 
ficients, r and  t,  above  can  represent  either  the  parallel  or  the  per- 
pendicular components.  The  subscripts  on  r and  t describe  the  boundary 
the  light  beam  hits,  for  example,  ti0,  is  the  transmission  coefficient 
for  light  going  from  medium  1 to  medium  0. 

The  total  amplitude  for  all  reflected  light  is 


R = Rj  + R2  + R3  + R4  +. 


(2.52) 


Using  the  pattern  developed  in  Eqs.  (2.49)  and  (2.50),  one  can  write 


R - r + 

S'  " r°i  + r 


10  1=  1 


r-  r -2  i 5 -1  n 

Lr12rloe  ] 


(2.53) 


This  sum  is  just  a geometric  series  so  that  it  is  trivial  to  obtain  a 
closed  form.  Also,  certain  simplifications  are  possible.  From  the 
explicit  expressions  shown  in  Eq.  (2.28),  one  can  deduce  that 


1 0 


= -r 


0 1 


and 


t t 
0110 


0-r,\> 


(2.54) 


» 


rr 
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Using  Eqs.  (2.54)  and  performing  the  sum  in  Eq.  (2.53),  one  obtains: 


R = roi  + riz  e"£1° 

A 1 + r01r12  e_2i6 


(2.55) 


For  the  separate  parallel  and  perpendicular  components  of  the 
light,  one  has: 


Rc  „s  s -i 26 

„ _ _s  _ r0 1 + rt2  e 

s Ac  i „s  -i25 

s 1 + r10  ri2  e 


. RP  rg,  + r?,  e-i2i 

" #P  1 ♦ r»,.r?2  e-12{ 


(2.56) 


so  that  Eq.  (2.47)  becomes: 


tan  e 


+iA  _ rn i + r?7  e~l2lS  1 + r?„  rf2  e~l26 
1 + r?„  r?2e1'-26  r!0  + r?2  e'i26 


(2.57) 


Eq.  (2.57)  is  the  desired  connection  between  the  parameters  and 
A,  determined  with  the  ell ipsometer,  and  the  media  parameters 


n0  e 1 = index  of  refraction  of  air. 


hf  = nf  (1  + iKf), 


(2.58) 


»,  ="s  (1  + 1Ks>- 


d = film  thickness 


characterizing  the  film  and  substrate.  Since  these  media  parameters 
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enter  Eq.  (2.57)  through  expressions  of  the  form  given  in  Eq.  (2.28), 
it  is  clear  that  the  relationship  between  ^ and  A and  the  media  para- 
meters is  arithmetically  rather  complicated.  Nevertheless,  if  only  two 
of  the  media  parameters  are  unknown,  then  measurement  of  ip  and  A for  a 
given  angle,  e^,  permits  computation  of  the  two  unknowns.  This 
computation  is  accomplished  by  digital  computer. 

There  exists  a vast  collection  of  literature  on  the  subject  of 
ellipsometry.  The  Bibliography  attached  presents  a useful  sampling  of 
this  previous  work.  This  literature  deals  not  only  with  the  relation- 
ship of  the  index  parameters  ns,  K$,  n^,  and  t to  the  quantities 
and  A,  but  also  with  the  alignment  and  use  of  the  instrument  itself. 

The  next  section  relates  the  values  of  ip  and  A to  the  actual  method  by 
which  these  parameters  are  obtained. 


SECTION  III 


THE  MECHANICS  OF  ELLIPSOMETRY 

The  ellipsometer  essentially  consists  of  two  polarizers  and  a 
quarter  wave  plate.  Light  is  supplied  by  a HeNe  laser  and  a power 
meter  is  used  as  a detector. 

The  light  from  the  laser  passes  through  a collimator  to  a polar- 
izer then  through  the  quarter  wave  plate  to  the  sample.  After 
reflecting  off  the  sample  the  light  wave  passes  through  the  second 
polarizer  to  the  detector.  Fig.  (3.1)  gives  a schematic  representation 
of  the  ellipsometer  in  the  plane  of  incidence.  The  collimator  arm  with 
its  polarizer  and  quarter  wave  plate  is  stationary.  The  telescope  arm 
with  its  polarizer  can  be  rotated  about  a center  axis.  A table  on 
which  the  sample  is  to  be  placed  is  centered  at  this  axis. 

The  angle  of  the  telescope  arm  can  be  read  directly  from  a scale. 
This  angle  is  measured  from  the  collimator  arm  as  0°,  so  that,  when 
light  shines  straight  through  from  the  collimator  into  the  telescope 
arm,  as  in  Fig.  (3.2),  the  angle  is  180°. 

The  angle  of  incidence,  e-,  can  be  found  using  this  angle  reading 
for  a sample  centered  on  the  table.  The  sample  table  can  be  tilted  to 
adjust  the  sample  so  the  plane  of  incidence  and  the  plane  formed  from 
the  two  arms  of  the  ellipsometer  are  the  same. 

The  polarizers  on  both  the  collimator  and  telescope  arms  are 
mounted  on  a circle  calibrated  to  .01°.  They  are  set  to  read  0°  when 
the  electric  vector  of  the  light  beam  is  horizontal  (in  the  plane  of 
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Fig.  3.2  Schematic  Representation  of  the  Ell ipsometer  with 
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incidence).  A quarter  wave  plate  is  essentially  a flat  transparent 
window  with  two  axes,  in  the  plane  of  the  window  at  approximately  90° 
to  each  other,  called  the  fast  axis  and  the  slow  axis.  For  light 
passing  through  it,  the  quarter  wave  plate  provides  a phase  delay  of 
90°  for  the  portion  of  the  light  polarized  along  the  slow  axis  with 
respect  to  the  portion  of  the  light  polarized  along  the  fast  axis.  The 
direction  of  polarization  of  light  is  considered  to  be  the  direction  of 
the  electric  field  vector. 

The  ellipsometer  is  used  by  adjusting  the  two  polarizers  so  that 
light  to  the  detector  is  a minimum.  This  is  done  by  adjusting  the  com- 
bination of  the  polarizer  and  quarter  wave  plate  in  the  collimating  arm 
so  that  a delay  between  the  parallel  and  perpendicular  components  is 
created  which  is  just  opposite  that  caused  from  the  reflection  of  light 
from  the  sample.  The  light  reflected  from  the  sample  is  then  linearly 
polarized  and  by  placing  the  second  polarizer  at  an  angle  of  90°  the 
light  can  be  extinguished. 

It  is  useful  to  determine  the  relationship  between  the  amplitudes 
A and  A , defined  in  Section  II,  and  the  angle  P of  the  polarizer.  The 

r •* 

quarter  wave  plate  is  set  with  the  fast  axis  at  an  angle  of  a = 45° 
with  respect  to  the  (horizontal)  ray  plane.  To  be  specific,  a unit 
vector  along  the  fast  axis  and  one  along  the  slow  axis  are  defined: 


6fast^[S  + e(1)] 

6slo»s^t5-P(i)] 


(3.1) 


fe 
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The  polarizer  may  be  set  at  some  arbitrary  angle  P with  respect  to  the 
ray  plane.  A unit  vector  identifying  the  direction  of  polarization  of 
light  from  the  polarizer,  can  be  defined: 


epoi  = (cosP)  p^  + (sinP)  £ 


(3.2) 


The  input  light  comes  out  linearly  polarized  at  angle  P (polarized 
parallel  to  epQ^ ) when  it  passes  through  the  polarizer.  If  the  magni- 
tude of  the  amplitude  of  the  electric  vector  is  E0,  then  the  parallel 
and  perpendicular  components,  with  respect  to  the  ray  plane,  after 
passing  through  the  polarizer  are  (in  the  notation  of  Section  II): 


Ep  = E0  cosP  e 


Es  = Efl  sinP  e 


(3.3) 


The  light  then  passes  through  the  quarter  wave  plate  and  the  components 
along  the  fast  and  slow  axes  become: 


Efast  = E0[cosP  cos45°  + sinP  sin45°]  e~lwt 
Eslow  = Eo^cosP  cos135°  + sinP  cos45°]  e_1 


(3.4) 


where  8 is  the  phase  delay  from  the  slow  axis  with  respect  to  the  fast 
axis.  From  these  components  the  parallel  and  perpendicular  components 


can  be  projected: 


£(i)  = A s + A p(i) 
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where 


A$  = ^-  [(cosP  + sinP)  + (sinP  - cosP)  ei8]e"iu)t 
Ap  = ^-  [(cosP  + sinP)  - (sinP  - cosP)  e^8]e*’(i)t 


(3.5) 


If  the  wave  plate  is  adjusted  to  be  a quarter  wave  plate  for  the 
wavelength  of  the  laser  source,  then  the  delay  B is  90°  and  these 
equations  can  be  simplified  to: 


A 


s 


A 


P 


-i(u)t  - P + 45°) 
-i  (tot  + P - 45°) 


(3.6) 


Eqs.  (3.6)  give  the  desired  relationship  between  A$  and  Ap  and  the 
angle  P of  the  polarizer. 

The  angles  \ Ij  and  A introduced  in  Eq.  (2.47)  are  defined: 

. A M 

(taml))elA  = -fjf-  . (3.7) 

s p 

In  order  to  experimental ly  determine  ij;  and  A,  the  following  strategy 
can  be  employed.  One  can  vary  the  angle  P of  the  polarizer  until 
linearly  polarized  light  is  reflected  from  the  sample.  This  condition 
can  be  detected  when  there  exists  an  angle  A of  the  analyzer  (the 
second  polarizer)  which  will  completely  block  the  light  reflected  from 
the  sample. 

Thus,  when  the  detector  is  "nulled"  to  the  lowest  attainable 
reading,  then  the  light  reflected  from  the  sample  is  linearly  polarized. 


> 


Under  these  conditions  there  are  only  two  possibilities  for  the  phase 


of  Rp/R5: 


f RrJl 

(a)  Phase  SrT3  0 ’ in  whl'ch  case  one  can  take 

R A [ A 1 

tamp  = f-  and  A = ^ Phase^  £-> 

fR  1 

(T)  Phase  >=  ±180°»  in  which  case  one  can  take 


(3.8) 


tan>p  = ■ r and  A = -180°  + Phase<  ^ 


(3.9) 


The  expressions  for  tamp  shown  above  are  permitted  because  |AS/Ap|  = 1 
as  shown  in  Eq.  (3.6).  Also  from  Eq.  (3.6),  ones  finds 

Phase!  = -2P  + 90°  . (3.10) 


Phase<  = -2P  + 90c 


Using  Eq.  (3.10)  in  Eqs.  (3.8)  and  (3.9)  the  two  cases  may  be 


expressed 


R R 

(a)  ^ 0,  tamp  = ^ , and  A = -2P  + 90° 


(3.11) 


R R 

(IT)  ^0,  tamp  = , and  A = -2P  - 90°  . (3.12) 

Since  tne  angle  A of  the  analyzer  is  set  at  right  angles  to  the  plane 
of  polarization  of  the  light  reflected  from  the  sample,  it  is  easy  to 
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find  the  connection  between  A and  \p  as  illustrated  in  Fig.  (3.3): 


Rp  =V/Rp  + R*  cos (A  + 90°) 


Rs  * yRp  + R*  sin(A  + 90°) 


so  that,  from  either  Eq.  (3.11)  or  Eq.  (3.12), 


(3.13) 


cos(A  + 90°) 
sin  AT  90°) 


1=1-  ~ r 1=1-  tanA  | (3.14) 


Furthermore, 


% _ . .1  . J Case  @ pertains  if  tanA  < 0 

n - “Lana  > \ r /-m ^ r\ 


Case  (B)  pertains  if  tanA  > 0 . (3.15) 


Solutions  of  Eq.  (3.14)  for  \p  can  be  restricted  so  that 


ip  e [o,90°] 


(3.16) 


with  the  following  correspondences: 


A € [0,90°]  =>Case  (Jl)  and  = A 
A e [90°, 180°]  =*>  Case  (a)  and  <p  = tt  - A 

A e [-90°, 0]  =>  Case  (/\)  ard  \p  = -A 

A £ [-180°, -90°]  =?>  Case  (IT)  and  ip  = n + A 


(3.17) 


Without  restriction,  multiples  of  360°  can  be  added  to  or  subtracted 
from  the  measured  angle  A to  bring  it  into  one  of  the  ranges  considered 
in  Eqs.  (3-17).  Similarly,  it  is  clear  from  Eq.  (3.7)  that  A need  only 
be  determined  mod  360°.  Thus  multiples  of  360°  can  be  added  to  and 
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subtracted  from  the  angles  A determined  in  Eqs.  (3.11)  and  (3.12). 

Thus,  the  angle  A may  be  taken  to  lie  in  the  range 

A e [0,360°]  . (3.18) 


The  two  cases  encountered  above  arise  because  extinction  by  the 
analyzer  can  be  achieved  for  two  different  angles  P and  P',  corres- 
ponding to 


[-2P  + 90°]  = A , Case  (a) 

and  (3.19) 

[-2P  + 90°]  = A + 7f  , Case  (?)  . 

It  is  important  to  note  that  when  a sample  is  placed  in  the 
ell ipsometer,  there  are  two  different  sets  of  angles  (Pi,  Ai)  and 
( P 2 , A2)  that  can  be  determined  by  the  measurement  procedure.  Either 
A 1 or  A2  will  have  a negative  tangent  (Case  (X))  and  the  other  will 
have  a positive  tangent  (Case  (?)).  The  test  shown  in  Eq.  (3.15)  is 
used  to  identify  which  angle  set  is  Case  (X)  and  which  is  Case  (?). 

The  angles  4>  and  A obtained  from  Eq.  (3.14)  and  the  appropriate  choice 
of  Eq.  (3.11)  or  Eq.  (3.12)  are  exactly  the  same 


V*® 

4®“A© 


in  the  absence  of  experimental  errors.  The  use  of  both  angle  sets 
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(Pi,  Ai)  and(P2»  A2)  allows,  therefore,  an  experimental  cross  check. 
Both  sets  were  measured  and  used  in  the  experimental  work  described  in 
the  following  sections.  In  the  next  section  an  error  analysis  and 
sensitivity  study  for  ellipsometric  data  are  described. 


SECTION  IV 


AN  ERROR  ANALYSIS  AND  SENSITIVITY  STUDY  FOR  ELLIPSOMETRIC  PARAMETERS 

When  using  the  ellipsometric  measurements  of  and  A to  examine  a 
film  or  substrate,  it  is  important  to  have  some  idea  of  how  much  error 
and  A contain.  Also,  since  both  the  real  and  imaginary  part  of  the 
refractive  index  for  both  the  substrate  and  film,  as  well  as  the  film 
thickness  (denoted  by  t in  this  and  the  following  sections)  are 
dependent  upon  \p  and  A,  it  is  useful  to  know  how  much  error  in  these 
quantities  will  result  from  errors  in  \p  and  A. 

To  examine  the  errors  in  ip  and  A,  the  equipment  used  in  making  the 
measurements  must  be  examined.  A Gaertner  L-119  ellipsometer  was  used 
for  the  measurements  made  for  this  study.  This  instrument  has  two 
arms,  the  collimator  and  the  telescope,  both  with  verticle  circles 
holding  polarizers  which  can  be  read  to  an  accuracy  of  0.01°  and  with  a 
clamp  and  tangent  screw  for  setting.  Mounted  on  a third  circle  is  a 
Soleil  Babinet  Compensator  which  can  be  adjusted  for  use  as  a quarter 
wave  plate.  Centered  between  the  two  arms  of  the  ellipsometer  is  a 
sample  table  which  can  be  rotated  and  adjusted  in  height.  The  colli- 
mator arm  is  fixed  and  the  telescope  arm  rotates  on  the  same  axis  as 
the  sample  table.  The  two  arms  have  height  and  level  adjustments  for 
alignment.  The  angle,  29..,  between  the  two  arms  can  be  read  from  a 
dial  to  20  seconds  of  arc. 

A prototype  detector  from  Newport  Research  Corporation  was  used  in 
making  the  measurements.  This  instrument  is  sensitive  to  light  inten- 
sities ranging  from  10"7  to  10"1  watts.  The  light  source  was  a He-Ne 
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laser  operating  at  a wavelength  of  6328A.  This  laser  was  made  by 
C.  W.  Radiation,  model  S-203H,  and  has  a maximum  output  of  10~3  watts. 

The  alignment  of  the  arms  was  made  so  that  light  through  the  col- 
limator could  pass  straight  through  (no  reflection)  into  the  telescope 
arm  and  out  a pinhole  exit.  A high  quality  mirror  surface  was  then 
used  to  insure  that  the  alignment  was  also  correct  for  reflection  into 
the  telescope  arm  for  all  values  of  6^,  the  mirror  surface  first  being 
tested  to  insure  that  it  was  perpendicular  by  reflecting  the  light  beam 
back  on  itself  from  both  sides  of  the  mirror. 

The  alignment  of  the  polarizers  was  done  next.  When  light  shines 
at  the  Brewster  angle  on  a piece  of  optically  flat  glass  of  index  n, 
the  reflected  light  is  perpendicular  to  the  plane  of  incidence.  At 
this  angle,  described  by  0.  = eD: 

I D 

tanGg  = n , (4.1 ) 

the  reflected  and  transmitted  waves  are  at  right  angles  to  each  other 
and  the  reflected  wave  does  not  receive  any  energy  for  oscillations  in 
the  plane  of  incidence.  By  using  this  property  of  the  Brewster  angle, 
the  analyzer,  fixed  on  the  telescope  arm,  will  extinguish  the  light 
when  it  is  oriented  horizontally  (A  = 0).  The  alignment  of  the 
analyzer  can  be  accomplished  by  adjusting  the  calibrated  circle  to  read 
A = 0.00  when  the  light  is  extinguished.  The  polarizer  on  the  colli- 
mator arm  can  now  be  adjusted  using  the  aligned  analyzer  to  make 
adjustments.  For  this  adjustment  the  arms  of  the  ellipsometer  are 
placed  so  the  light  shines  straight  through  and  there  is  no  reflection. 
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The  polarizer  angle,  P,  is  set  equal  to  0.00°  when  the  analyzer  angle, 
A,  is  set  to  90.0°  and  the  light  is  extinguished. 

Once  the  angles  of  the  analyzer  and  polarizer  are  set,  the  amount 
of  error  in  the  alignment  can  be  determined.  First,  without  any  sample 
to  reflect  the  light,  the  angle  A is  set  and  the  angle  P found  which 
extinguishes  the  light  to  the  detector.  The  error,  e,  between  the 
setting  of  the  two  polarizers  can  be  defined  by  the  relation: 

P = A + 90°  + e . (4.2) 

For  the  alignment  of  the  ellipsometer  used,  this  error  was  determined 
to  be  e = 0.03°.  This  relation  is  plotted  as  Curve  1 in  Fig.  (4.1). 

A mirror  surface  was  then  used  to  reflect  the  light  to  the  analy- 
zer. The  angle  A was  set  at  various  angles  in  the  neighborhood  of 
0.0°  and  for  each  value  of  A the  angle  P was  found  which  extinguished 
the  light.  For  P equal  to  exactly  90°  and  for  A equal  to  exactly  0°, 
the  error  due  to  any  phase  change  upon  reflection  would  be  zero  because 
there  would  be  no  horizontal  component.  At  angles  close  to  these 
values  (±0.2°)  the  amplitude  of  the  horizontal  component  is  still  very 
small  (0.35%)  compared  to  the  vertical  component  and,  therefore,  any 
phase  change  would  be  negligible.  The  data  from  these  measurements  is 
also  plotted  on  Fig.  (4.1).  Curve  2 on  this  figure  is  a straight  line 
lease  squares  fit  to  this  data.  The  intersection  of  the  two  curves 
locates  the  readings  of  A and  P which  actually  correspond  to  0.0°  and 
90.0°,  respectively.  The  corrections  for  A and  P so  determined  were: 
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Ac  = A + AA,  AA  = -0.07° 
Pc  = P + AP,  AP  = -.04° 


(4.3) 


The  corrected  angles  are  used  in  calculating  i|>  and  A and  all  further 
references  to  A and  P will  refer  to  the  corrected  angles,  Ac  and  P . 

After  the  polarizers  were  both  aligned,  the  quarter  wave  plate  was 
inserted  and  aligned.  The  quarter  wave  plate  used  is  a Soliel  Babinet 
Compensator.  This  compensator  can  be  varied  in  thickness  so  it  can  be 
used  as  a quarter  wave  plate  over  a wide  range  of  frequencies.  The 
quartz  used  as  the  retarder  in  the  wave  plate  is  actually  in  two 
pieces.  The  outer  edges  are  parallel  but  the  cut  between  the  two 
pieces  is  diagonal.  The  thickness  of  the  quartz  can  be  varied  by 
sliding  the  two  pieces  of  quartz  along  each  other.  Fig.  (4.2a)  shows 
the  alignment  of  the  quartz  pieces.  To  adjust  the  wave  plate  in  the 
ellipsometer,  the  fast  axis  of  the  compensator  was  first  set  to  an 
angle  of  a = 45.0°.  To  do  this,  the  analyzer  was  set  to  45.00°  and  the 
polarizer  was  set  to  135.00°.  With  the  fast  axis  of  the  compensator 
at  an  angle  of  approximately  45°,  the  angle  was  adjusted  for  the  light 
to  the  detector  to  be  at  a minimum.  At  this  point  the  fast  axis  scale 
was  set  to  read  45.00°. 

With  the  angle  scale  set  for  the  quarter  wave  plate,  the  thickness 
or  phase  delay,  6,  was  tuned  for  a delay  of  90°.  Leaving  the  fast  axis 
at  45°,  the  polarizer  was  set  at  0.00°  and  the  analyzer  at  90.00°,  the 
light  to  the  detector  was  then  minimized  by  changing  the  thickness.  The 
thickness  of  the  quartz  was  varied  over  a wide  range,  giving  six 
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Fig.  4.2  The  Two  Pieces  of  the  Compensator  Showing 

a)  the  Direction  of  Movement  for  Adjustment 
of  Thickness 

b)  the  Alignment  Error  in  the  Axis 
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readings  for  which  the  light  was  at  a minimum.  At  these  points  the 
amounts  of  delay  due  to  the  thickness  is  an  odd  multiple  of  180°.  The 
dial  which  adjusts  the  quartz  thickness  is  calibrated  in  dimensionless 
steps  of  0.01  and  the  spacing  between  the  minimum  readings  is  6.24. 
Assuming  the  dial  to  be  linear,  the  reading  for  a delay  of  ±90°  could 
be  calculated  by  adding  1.56  to  any  of  the  readings  taken.  An  error 
of  0.07  on  the  dial  would  cause  an  error  of  1.0°  in  the  delay. 

The  accuracy  of  the  alignment  of  the  compensator  can  be  checked 
in  several  ways.  With  the  thickness  set  for  a 90°  delay  and  the  fast 
axis  set  at  45°,  if  the  polarizer  is  set  at  0.00°  the  light  passing 
through  the  analyzer  to  the  detector  should  be  at  the  same  intensity 
for  all  settings  of  the  angle  A.  In  this  case  the  light  leaving,  the 
wave  plate  is  circularly  polarized.  The  alignment  of  the  fast  axis 
can  also  be  checked  by  setting  both  the  polarizer  and  the  compensator 
to 45.00°,  setting  the  analyzer  to  135.00°  and  changing  the  thickness  of 
the  quartz.  If  all  three  elements  are  aligned  there  should  be  no 
change  in  the  nulled  intensity  of  the  light  to  the  detector.  For  this 
particular  compensator,  however,  the  two  pieces  of  quartz  were  slightly 
misaligned.  The  fast  axis  of  the  two  pieces  formed  an  angle  of  approxi- 
mately 0.20°,  see  Fig.  (4.2b).  This  could  be  shown  by  leaving  the  com- 
pensator at  45.0°,  setting  the  polarizer  at  135.00°,  and  adjusting  the 
analyzer  to  minimize  the  light  to  the  detector.  The  angle  A was  at 
45.09°  and  the  amount  of  error  was  +.09°  from  the  theoretical  angle  of 
45.00°.  Angle  A was  then  reset  to  45.00°  and  the  angle P was  adjusted  so 
that  a minimum  intensity  was  detected.  The  angle  measured  was  134.89°  so 
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that  the  amount  of  error  was  -0.11°.  The  conclusion  was  that  the 
first  piece  of  the  compensator  is  at  an  angle  of  44.89°  and  the  second 
piece  is  at  an  angle  of  45.09°  giving  the  difference  between  the  two 
pieces  as  0.20°. 

To  find  the  amount  of  error  this  will  cause  in  the  calculations  of 
ip  and  A,  some  theoretical  calculations  were  made.  A computer  program 
written  by  F.  L.  McCracken[5]  from  the  National  Bureau  of  Standards  was 
used  in  calculations  for  error  analysis  as  well  as  in  determining  ip  and 
A from  measurements  of  the  angles  A and  P.  Theoretical  values  of  A and 
P were  determined  for  a film  index  of  (nf  = 1.766,  K ^ = 0)  and  a sub- 
strate index  of  (n$  = 1.21  and  Ks  = 5.72),  and  film  thickness  varying 

o o 

from  0 to  600A  in  increments  of  100A.  Using  these  values  for  A and  P, 
the  values  of  \p  and  A were  calculated  for  conpensator  angles,  a,  of 
44.95,  45.0  and  45.05.  To  check  the  errors  in  \|>  and  A due  to  an  error 
in  the  delay  of  the  wave  plate  of  0.5°,  calculations  were  also  done 
using  a compensator  angle  of  45.0°  but  using  the  values  of  89.5°, 

90.0°,  and  90.5°  for  8.  The  results  of  these  calculations  are  shown  in 
Tables  (4.1)  through  (4.4).  Since  two  sets,  (Pi,  A])  and  (P2,  A2),  can 
be  found  for  each  ij>  and  A,  both  sets  were  used.  The  averages  of  the 
two  values  obtained  for  and  A were  also  calculated  and  denoted  by  the 
subscript  'avg'.  A computer  sensitivity  study  showed  that  ij>avg,  corres- 
ponding to  errors  in  a by  either  +0.05°  or  -0.05°,  agreed  with  the 
value  of  ip  obtained  with  no  error  in  a.  These  results  are  shown  in 
Table  (4.1)  and  illustrate  the  value  of  using  >J'avg-  for  an  error  in  8 
of  ±0.5°  the  average  was  only  a little  worse  giving  an  error  in  the 
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third  decimal  place  as  shown  in  Table  (4.2). 

For  A,  however,  the  errors  were  more  pronounced.  An  error  in  a 
of  ±0.05°  gives  an  error  of  ±0.1°  for  AaVg,  and  an  error  of  ±.5°  in  8 
gives  an  error  of  up  to  .006°  for  Aavg.  These  results  are  illustrated 
in  Tables  (4.3)  and  (4.4). 

The  slight  misalignment  of  the  quartz  pieces  in  the  compensator, 
mentioned  above,  produced  the  ambiguity:  a = 44.89°,  but  also 
a = 45.09°.  The  angle  44.89°  was  used  in  the  computer  code  for 
analysis  of  the  data.  Additional  checks  showed  that  the  0.20  discre- 
pancy would  produce  negligible  errors  in  \|>  . 

The  errors  produced  in  x p and  A by  errors  in  A and  P can  be 
expressed  simply  since 

= A ‘ or  it  + A' 

for  Case  (b) 

A = -2P1  - 90° 

and 


xp  - -A  or  it  - A 
A = -2P  + 90° 


for  Case 


If  the  amount  of  error  in  the  angle  A is  e,  then  the  error  in  i|>  would 
be  e if  a single  value  of  A were  used.  But  if  two  values  of  A were 
found,  say,  one  corresponding  to  A between  0°  and  90°,  and  the  other 
for  A between  -90°  and  0°  then  the  error  is  reduced  because  some 
sources  of  error  give  cancelling  effects.  For  example,  an  error  in  the 
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zero  position  of  the  scale  for  A will  cancel 

- (A1  +_  e)  + 1.  -A  + e 1 _ 

^avg  2 M 

Unfortuantely , an  error  in  the  zero  position  of  the  scale  for  P does 

not  average  out  for  errors  in  A,  . If  the  error  in  P is  6,  the 

dvg 

amount  of  error  in  A can  be  calculated  by 

A _ 2 (P * + 6)  - 90°}  + {-  2 (P  + 6)  + 90°} 

aavg  “ 2 

A A-  26  + A'  - 26  . „r 

4avg  * 2 ' 4 - 25  ■ 

The  error  in  A and  in  Aavg  is  double  that  in  the  measurement  of  P or  P1. 

Another  and  more  critical  source  of  error  arises  from  the  method 
in  which  the  sample  is  aligned  in  the  ellipsometer.  Although  the  cal- 
culations determining  the  values  of  iJj  and  A from  angles  A and  P do  not 
depend  on  the  angle  of  incidence,  the  use  of  and  A in  determining  the 
values  of  the  film  thickness  or  the  index  of  refraction  is  dependant  on 
this  angle.  A small  angle  of  incidence  error  can  cause  large  errors  in 
ip  and  A.  The  placement  of  the  sample  is  critical  because  an  error  of 
5mm  from  the  center  of  the  table  can  create  an  error  of  0.5°  in  the 
angle  of  incidence,  0...  Fig.  (4.3)  shows  the  values  of  and  A for  a 

o o 

film  with  thickness  varying  from  0A  to  200A.  The  film  index  is  again 
taken  to  be  (nf  = 1.76,  K = 0)  and  the  substrate  index  is  n$  = 1.21 
(1  + j5 . 72 ) . 

Because  this  source  of  error  was  recognized  as  the  quantity  which 
would  most  likely  limit  the  resolution  of  the  ellipsometer,  a systematic 


procedure  for  aligning  each  sample  was  used.  Pinholes  on  each  end  of 
the  telescope  tube  were  carefully  aligned.  If  a sample  were  carefully 
placed  at  the  center  of  the  sample  table,  the  reading  on  the  scale  for 
(26..)  would  be  accurate  to  ±20  sec  of  arc.  The  light  would  then  pass 
through  the  center  of  the  pinholes.  For  a misaligned  sample,  the 
light  would  not  pass  straight  through  the  telescope  tube,  but  would  be 
inclined  by  an  amount  66  as  indicated  in  Fig.  (4.4).  The  maximum  angle 
of  error  can  be  defined  by  the  diameter  of  the  pinhole  (.2cm)  and  the 
length  of  the  telescope  tube  (28cm). 


66  = Tan'1  (5!^) 

max  '28cm' 


0.4° 


The  error  in  the  angle  of  incidence  is  half  this  value  giving  a possi- 
ble error  of  ±0.2°.  The  alignment  of  a sample  was  always  checked  in 
two  configurations  as  shown  in  Fig.  (4.4).  This  procedure  served  as  a 
partial  double  check  on  the  pinhole  alignment,  also,  because  the  check 
in  two  configurations  would  be  impossible  if  the  pinholes  are  not  well 
al igned. 

An  added  benefit  from  using  this  method  to  align  the  sample,  is 
that  the  sample  tilt  is  also  checked.  The  sample  cannot  pass  this  test 
if  it  leans  forward  or  backward  more  than  0.2°.  This  tilt  can  cause  an 
error  in  \p  for  the  two  sets  of  measurements,  but  for  ij/  the  error  is 
negligible.  The  tilt  angle  also  causes  an  error  of  less  than  .01°  in 

Aavg‘ 

Another  source  of  error,  that  due  to  the  transmission  properties 
of  the  wave  plate,  has  also  been  examined.  A transmission  difference 


of  5%  between  the  fast  and  slow  axes  would  cause  large  discrepancies 
in  A (about  5°)  between  the  two  sets  of  measurements  at  each  point. 
The  average  of  these  two  measurements  gives  a very  small  error  of 
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0.03°  in  A,  . The  error  in  to,  „ is  also  about  0.03°.  Since  the 

avg  avy 

maximum  discrepancy  in  the  sets  of  measurements  actually  taken  was 
0.66°  for  Aau_  and  .57°  for  \ any  error  due  to  the  transmission 

avy  avg 

properties  of  the  wave  plate  would  have  to  be  negligible. 

The  conclusion  drawn  from  this  study  of  errors  is  that  the  error 
in  the  angle  of  incidence  has  the  largest  effect.  The  angle  of  inci- 
dence is  correct  to  the  accuracy  of  ±0.2°  and  the  effect  of  this  error 
on  values  of  ^ or  A must  be  examined  for  the  particular  set  of  values 
given  for  the  film  and  substrate  properties  to  determine  expected 
error  bounds. 

The  precision  of  the  measurements  can  be  estimated  by  examining 

the  data.  Two  sets  of  measurements,  (Pi,  Ax)  and  (P2,  A2),  were  used 

to  obtain  the  values  of  A.w_  and  to.  At  each  point  on  the  sample, 

avg  avg 

the  numbers  A*,  A2,  and  AflVg  were  determined.  The  maximum  discrepancy 
in  A was  0.66°  and  the  average  discrepancy  was  0.21°.  Thus  the  average 
precision  of  the  procedure  for  A was  ±0.11°  to  ±0.33°.  The  maximum 
discrepancy  in  was  0.57°  and  the  average  discrepancy  was  0.18°.  Thus 
the  average  precision  of  the  procedure  for  ijj  was  ±0.09°  to  ±0.29°. 

Once  some  idea  is  formed  concerning  the  amount  of  error  in  \p  and 
A,  it  is  necessary  to  relate  this  error  to  the  error  produced  in  the 
index  of  refraction  and  in  the  film  thickness.  This  is  done  by  use  of 
curves  of  constant  index  and  thickness  plotted  against  the  values  of  \J> 
and  A. 
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The  simplest  case  to  study  is  that  of  a bare  substrate,  since 
there  are  only  two  media  parameters  determined  by  ip  and  A.  These  two 

variables  are  the  real  part,  ns  and  the  imaginary  part,  K$,  of  the 

index  of  refraction,  where 

n = ns  (1  + ifC$)  . 

Fig.  (4.5)  shows  curves  of  constant  n$  and  K$  on  a and  A plot.  It 
should  be  noted  that  for  <s  = 0,  there  is  no  phase  change  upon  reflec- 
tion and  A has  a value  of  0.  As  the  imaginary  part  of  the  index 
increases  the  phase  between  the  reflected  parallel  and  perpendicular 
components  increases  and,  thus,  the  value  of  A increases. 

Fig.  (4.6)  shows  the  theoretical  values  of  ^ and  A for  a film  of 

nf  = 1.6,  Kf  = 0 on  a substrate  of  n$  = 1.21  with  no  absorption.  The 

° O 

thickness  varies  from  0 to  2700A.  At  a thickness  of  238 1A  the  curve 
closes  and  starte  repeating,  this  happens  because  the  film  thickness 
produces  in  a phase  shift  of  360°  at  this  point.  For  a film  with  no 
absorption  the  amplitude  and  phase  would  be  the  same  at  the  surface  for 
a phase  shift  of  any  number  of  complete  cycles.  The  thickness  of  any 
film  to  cause  a full  360°  phase  shift  can  be  calculated  as: 

t = Xq _ \ 

c 2n  cos  (arcsin  (sinG^/n)}  2 /n^  - sin*©- 

This  thickness  is  a function  of  the  wavelength  of  the  light  A0,  the 
film  index  n,  and  the  angle  of  incidence  0... 

Fig.  (4.7)  shows  a family  of  curves  generated  by  varying  the  film 
thickness,  t.  These  curves  are  for  the  film  index,  n^,  varying  from 
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1.1  to  1.8.  The  substrate  is  the  same  as  in  the  previous  figure.  It 

can  be  noted  that  as  the  film  index  approaches  the  substrate  index  the 

curve  formed  becomes  smaller  and  at  the  point  where  n^  » ng  only  a 
single  point  is  found.  This  point  can  be  compared  to  a point  with  the 

same  ij;  and  A in  Fig.  (4.5),  and  it  can  be  noted  that  this  is  the  value 

for  a bare  substrate  with  that  index. 

Using  a film  index  of  n^  = 1.766  and  the  real  part  of  the  sub- 
strate index  to  be  n&  = 1.21,  the  imaginary  part  of  the  substrate  Ks 
was  increased  from  0.0  to  5.72.  Fig.  (4.8)  shows  thickness  curves 
corresponding  to  this  range  of  K$.  All  of  the  thickness  curves  are 
cyclic.  The  zero  thickness  trajectory  as  a function  of  Ks  is  traced 
by  the  dotted  curve. 

Since  the  material  used  in  this  study  has  a substrate  of  aluminum, 
Fig.  (4.9)  shows  a family  of  curves  with  the  film  index  changing  from 
1.1  to  1.8  on  a substrate  Of  n$  = 1.21  and  K$  = 5.72,  values  similar  to 
aluminum. 

Fig.  (4.10)  shows  the  change  in  the  curves.  As  the  real  part  of 
the  substrate  index,  ns,  is  varied  from  1.0  to  2.0.  The  imaginary 
part  of  the  substrate  index,  Ks  = 5.72,  and  the  film  index,  n^  = 1.766, 
are  constant.  Using  Figs.  (4.8)  and  (4.10),  it  can  be  seen  that  a 
small  error  in  the  values  assumed  for  real  and  imaginary  part  of  the 
substrate  index  will  not  cause  a large  shift  in  the  values  determined 
for  the  film  index  and  thickness. 

When  the  film  over  the  substrate  becomes  absorptive,  the  curves 
are  no  longer  cyclic.  Fig.  (4.11)  shows  the  curves  formed  by  varying 
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Fig.  4.10  Variation  in  Thickness  Curves  for  a Change  in  n 


Fig.  4.11  Variation  in  Thickness  Curves  as 
Kf  Varies  from  0 to  .2 


the  thickness  of  the  film  on  a substrate  of  aluminum.  The  real  part 
of  the  film,  nf  = 1.76,  is  constant,  but  the  imaginary  part,  Kf,  is 
considered  for  the  three  values:  0.0,  0.02,  and  0.2.  The  curve  formed 
with  Kf  = 0.2  is  shown  as  the  thickness  gets  large.  As  this  happens, 
enough  light  is  absorbed  that  the  reflected  light  sees  only  the  film 
and  the  values  of  ip  and  A approach  those  for  a substrate  with  the  index 
of  the  film.  This  point  can  be  compared  with  the  substrate  value  in 
the  enlarged  view  in  Fig.  (4.12). 

Since  the  curves  formed  from  a smaller  imaginary  value,  = 0.02, 
do  not  change  as  fast  as  for  = 0.2,  they  were  not  plotted  for  larger 
values  of  the  thickness.  At  a much  thicker  film,  these  curves  also 
approach  the  ip  and  A values  of  a substrate  with  the  index  of  the  film. 

Now  that  the  various  sources  of  error  and  parameter  sensitivities 
have  been  explored,  the  general  description  of  the  research  tool  used 
in  this  study  is  complete.  In  the  next  section,  the  use  of  ellipso- 
metry  for  the  study  of  oxide  films  on  pure  aluminum  will  be  described. 
Because  such  films  have  certain  unpleasant  properties,  certain 
strategies  must  be  developed  in  order  to  gain  useful  results  from 
ellipsometry.  The  strategies  of  the  method  are  generalized  even  more 
in  Section  VI  in  an  attempt  to  use  ellipsometry  to  study  oxide  layers 
on  Alclad. 
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Fig.  4.12  Enlarged  View  for  the  Changes  in  n$  and  K$ 
for  a Bare  Substrate 
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A STUDY  OF  OXIDE  FILMS  ON  PURE  ALUMINUM 

In  order  to  attempt  to  develop  a method  for  measuring  the  thick- 
ness of  oxide  films  on  the  rough  surfaces  of  alclad  materials,  oxide 
films  were  studied  first  on  optically  flat  surfaces  of  pure  aluminum. 
Samples  of  optically  flat  fused  silica  covered  with  an  evaporated  film 
of  99.999%  pure  aluminum  were  used  for  this  study.  In  air  an  oxide 
film  grows  very  fast  on  a bare  aluminum  surface.  R.  W.  Fane  and 
W.  E.  J.  Neal  [6]  have  monitored  changes  of  A over  a period  of  100  hours. 
The  results  are  given  in  Fig.  5.1  and  show  a linear  relationship  be- 
tween the  change  in  A and  the  log  of  time.  A change  in  the  value  of  A 
by  100  minutes  corresponds  approximately  to  a 12$  change  in  thickness. 
They  also  determined  that  the  final  thickness  of  the  oxide  film  was  be- 
tween 40  and  55$  depending  on  the  initial  state  of  oxidation,  time  of 
exposure,  and  humidity. 

J.  H.  Halford,  F.  K.  Chin,  and  0.  E.  Norman[7]  state  that  this 
oxide  layer  grows  slowly  in  a vacuum,  but  accelerates  when  exposed  to 
air.  They  estimate  the  asymptotic  thickness  to  lie  between  40  and  50$. 
Exposure  of  the  oxide  layer  to  0 ^ or  moist  0£  accelerates  the  growth. 

M.  J.  Digman  [8]  studied  the  oxide  growth  at  the  elevated 
temperatures  of  250°,  350°  and  450°C.  The  films  were  heated  and  at 
intervals  the  barrier  voltage  was  measured.  The  barrier  voltage  is 
proportional  to  the  oxide  thickness.  Table  5.1  shows  the  results.  The 
lowest  temperature,  250°C,  produced  a small  change  over  a 600  hour 
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TABLE  5.1 

RELATIONSHIP  OF  HEATING  TEMPERATURE  AND  TIME 
AS  RELATED  TO  THE  BARRIER  VOLTAGE 


Oxidation 

Temp. 


♦ 


450°C 

* 


Time,  Hr 
(per  cent) 
0.0667 
0.25 
1.0 
4 
16 
24 


Barrier 

Voltage 


1.78 

2.20 

2.47 

3.29 

3.66 

2.89 


1.0  1.70 

4 2.04 

350°C  16  .234 

64  2.83 

256  3.36 


1.0  1.33 

5 1.62 

250°C  24  1.75 

120  2.02 

600  2.37 
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This  data  was  taken  from  Dignam[8] 


period,  while  for  the  higher  temperatures  the  change  over  smaller 
periods  was  greater.  At  450°C  after  24  hours  there  is  a sudden 
decrease  in  the  barrier  voltage.  This  is  due  to  a change  in  the  oxide 
structure.  Films  formed  on  aluminum  below  500°C  are  either  amorphous 
material  or  material  composed  of  extremely  small  crystals.  At  higher 
temperatures  gamma  alumina  crystallites  are  formed.  No  evidence  of  the 
crystallites  were  found  for  oxidation  periods  up  to  256  hours  at  350°C. 
Digman  also  studied  samples  anodically  oxidized  and  annealed  at  450°C 
and  found  them  to  have  the  same  conduction  properties  as  the  samples 
with  the  oxide  grown  in  the  250°  to  450°C  range.  He  concluded  that  the 
film  structures  were  apparently  the  same.  The  oxide  film  formed  on 
aluminum  is  transparent  and  has  a refractive  index  of  1.65. 

In  this  study  the  samples  were  heated  in  an  oven  to  350°C  and  kept 
at  this  temperature  for  a period  of  time  of,  typically,  two  hours.  The 
samples  were  allowed  to  cool  overnight  before  the  surface  was  examined 
with  an  ell ipsometer. 

The  first  sample  studied.  Sample  A,  had  a IOkA  layer  of  aluminum. 
Five  different  points  on  the  sample  were  marked  so  that  measurements 
could  be  made  at  the  same  spot  each  time.  These  five  spots  were  num- 
bered 1,  2,  3,  4,  5 for  reference  purposes.  The  sample  was  heated  16 
times  and  Table  5.2  shows  the  measurements  made  prior  to  the  first 
heating  and  after  each  successive  heating.  The  number  of  hours  shown 
in  Column  1 is  the  total  number  of  hours  at  350°C,  but  excludes  the 
heat  up  and  cool  down  time.  For  each  spot  on  the  sample  this  data  was 
plotted,  showing  that  the  consecutive  measurements  generate  a A vs.  \p 
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curve.  These  data,  shown  in  Figs.  (5.2)-  (5.6),  all  follow  approximate- 
ly the  same  trajectory. 

Another  sample.  Sample  B,  with  a 500&  film  of  aluminum  was  also 
tested.  This  sample  also  was  vacuum  deposited,  but  was  covered  with  a 
layer  of  collodion  to  prevent  oxidation  of  the  aluminum.  After  re- 
moving the  collodion,  measurements  were  taken  with  the  ell ipsometer. 

The  sample  was  then  heated  twice  with  measurements  being  taken  after 
each  heating  at  two  distinct  spots,  a and  b,  on  the  sample.  Fig.  (5.7) 
shows  these  data  in  relation  to  those  taken  at  Point  1 on  Sample  A 
(from  Fig.  (5.2)). 

A study  was  then  undertaken  in  an  attempt  to  relate  the  values  of 
1 1>  and  A on  these  curves  to  the  index  of  refraction  of  both  the  sub- 
strate and  film  and  to  the  film  thickness.  Fig.  (4.9)  in  Section  IV 
shows  curves  for  a substrate  index  in  the  range  of  aluminum  for 
varying  thickness  and  film  index.  Fig.  (5.8)  shows  an  enlarged  view 
for  the  same  substrate  index,  but  with  film  thickness  between  0 and 
35oX.  The  combined  data  for  the  five  spots  on  Sample  A and  the  two 
spots  on  Sample  B are  shown  on  the  figure.  The  change  in  the  film  in- 
dex between  1.6  and  1.7  has  a negligible  effect  on  the  curve  fit. 
Comparison  of  the  curve  to  the  measured  data  shows  that  the  measure- 
ments lie  farther  and  farther  to  the  left  as  the  thickness  increases. 

An  effort  was  made  to  fit  a curve  to  the  data  by  assuming  that  the  film 
had  an  effective  index  that  was  complex.  The  assignment  of  such  a com- 
plex index  can  be  motivated  as  follows.  Since  a rough  surface  scatters 
some  of  the  light,  the  light  detected  from  the  reflection  is  reduced 
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Fig.  5.2  The  Data  Measured  on  Sample  A at  Spot  1 
for  Successive  Heatings 


I 


130 


125 


120 


115 


110 


105 


100 


95 


41 


I 


Fig.  5.4  The  Data  Measu 
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Fig.  5.5  The  Data  Measured  on  Sample  A at  Spot  4 
for  Successive  Heatings 
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Fig.  5.8  Comparison  of  Data  to  Curves  for 
nf  = 1.6  and  1.7  (No  Extinction 
Coefficient) 
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by  the  loss.  The  extinction  coefficient,  Kf,  is  a term  accounting  for 
light  absorbed  in  the  material  and,  therefore,  not  reflected.  Thus, 
it  may  be  reasonable  to  treat  scattering  losses  effectively  as  absorp- 
tion losses.  Fig.  (5.9)  shows  a curve  superimposed  over  the  data, 
using  a film  index  of  fif  = 1.65  (1  + j0.06).  The  substrate  index  for 
this  curve  was  1.15  (1  + j5.75).  This  curve  fits  the  data  much  better, 
but  the  data  and  the  curve  tend  to  bow  in  opposite  directions.  An 
effort  to  make  a better  match  involved  making  the  imaginary  part  of 
the  film  index  directly  related  to  the  film  thickness.  Fig.  (5.10) 
shows  two  curves  using  the  model 

Kf  = (O.O! ) (5.1) 

where  t is  the  oxide  film  thickness  in  A.  The  substrate  indices  for 
these  two  curves  are  1.15  (1  + j5.45),  as  in  Fig.  (5.9),  for  one  and 
1.65(1  +J4.75)  for  the  other.  For  the  1.15  (1  + j 5 . 45 ) case  the 
standard  deviation  of  the  sample  data  was  evaluated  for  the  fixed 
value  Kf  = 0.06,  used  in  Fig.  (5.9),  and  for  the  Kf  model  of  Eq.  (5.1) 
shown  in  Fig.  (5.10).  The  deviation  was  3.39  for  the  curve  employing 
Kf  = 0.06  and  was  2.07  for  the  curve  based  on  the  Kf- thickness  model  of 
Eq.  (5.1).  The  arbitrary  unit  of  these  standard  deviations  was  1/20  of 
an  inch.  These  calculations  show  that  the  curve  in  which  the  absorp- 
tion term,  Ks,  varies  matches  the  data  much  better. 

C.  A.  Fenstermaker  and  F.  L.  McCrackin[9]  have  constructed 
theoretical  models  for  errors  in  ellipsometric  measurements  due  to 
surface  roughness.  They  modelled  surfaces  as  square  ridges,  triangular 
ridges,  and  pyramids  of  equal  height  and  width.  The  size  of  this 
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roughness  was  varied  from  0 to  SOOK.  The  model  was  applied  to  several 
materials  with  the  true,  complex,  indices  of  refraction  as  given  below: 


Material 

n 

K 

glass 

1.50 

0.0 

silicon 

4.05 

0.0068 

chrome 

3.00 

1.4 

mercury 

1.485 

3.061 

gold 

0.35 

7.0 

silver 

0.18 

19.0 

Since  the  true  index  of  glass  has  no  absorption,  it  is  assumed 
that  this  most  nearly  resembles  the  aluminum  oxide.  Therefore,  it  is 
these  results  which  would  be  interesting  for  this  case.  Fig.  (5.11) 
shows  the  calculated  changes  in  n and  K for  glass  due  to  surface  rough- 
ness obtained  by  Fenstermaker  and  McCrackin.  The  apparent  value  for 
the  real  part  of  the  index  changes  very  little,  but  the  apparent  value 
of  K increases  significantly  with  surface  roughness. 

Since  the  aluminum  oxide  layer  would  tend  to  grow  unevenly,  one 
might  expect  the  surface  to  become  rougher  as  the  thickness  increases. 
The  results  of  Fenstermaker  and  McCrackin  do  show  an  increase  in 
as  the  film  thickness  increases.  Fig.  (5.10)  confirms,  for  this  case, 
that  the  roughness  (as  measured  by  t)  appears  to  increase  the  effective 
value  of  the  extinction  coefficient,  K^. 

As  the  two  curves  in  Fig.  (5.10)  illustrate,  various  pairs  (ns, 

K ) can  be  used  to  achieve  a good  fit  to  the  data,  using  the 
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Fig.  5.11  Calculated  Changes  in  n and  K Due  to  Surface 
Roughness  of  Glass 


K-thickness  model  of  Eq.  (5.1).  Curves  resulting  from  other  suitable 
index  choices  are  shown  in  Fig.  (5.12).  This  lack  of  uniqueness  could 
be  resolved  if  the  actual  film  thickness  were  known.  As  shown  in 
Fig.  (5.12),  the  thickness  scales,  t,  on  the  curves  for  different 
(n$,  K$)  are  pairs  are  quite  different. 

Previous  determinations  of  the  substrate  refractive  index  for  pure 
aluminum  with  thin  films  show  a wide  range  of  variation. [6,7,10,1 1]  Some 
of  the  values  found  for  the  substrate  refractive  index  in  the  last 
twenty-five  years  are  shown  in  Fig.  (5.13).+  In  this  figure  the  real 
part  n$  is  shown  plotted  against  the  extinction  coefficient  Ks.  The 
curve  marked  A in  Fig.  (5.13)  gives  the  locus  of  index  values  which 
fit  the  data  of  the  present  study  when  the  fixed  value  K.p  = 0.06  is 
used,  as  in  Fig.  (5.9).  The  curve  marked  B in  Fig.  (5.13),  on  the 
other  hand,  gives  the  locus  of  index  values  which  fit  the  data  when  the 
model  of  Eq.  (5.1)  is  used,  as  in  Figs.  (5.10)  and  (5.12).  Curve  B 
is  specified  by  the  linear  relation 

K$  = 7.175  - |ns  . (5.2) 

It  may  be  seen  that  Curve  B provides  a nice  fit  to  the  widely 
scattered  determinations  by  other  workers.  That  Curve  B gives  a good 
fit  to  this  data  has  two  important  consequences.  First,  this  fit  gives 


+The  point  marked  F on  the  figure  is  only  an  estimate  made  for  6328& 
light,  based  on  the  work  of  Fane  and  Neal [6]  done  at  a wavelength  of 
5490&.  The  estimate  was  made  using  the  trend  of  the  data  with  A deter- 
mined by  Hass  and  Waylonis . [1 1 ] 


A - Halford,  Chin  and  Norman 
B - Schulz 

C - Hass  and  Waylonis 

D - Winterbottom 

E - Ehreneich 

F - Fane  and  Neal 

(Estimated  for  6328/1) 

G - Nyce  and  Skolnick 
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an  independent  indication  that  the  model  of  Eq.  (5.1)  may,  indeed  be 
very  suitable  for  the  oxide  films  on  pure  aluminum.  Secondly,  Curve  B 
offers  a way  of  reconciling  the  scattered  results  of  previous  workers. 
The  missing  piece  of  data,  at  least  in  the  present  study,  that  makes  it 
difficult  to  determine  which  point  on  Curve  B should  be  taken,  is  know- 
ledge of  the  actual  thickness  of  any  of  the  oxide  layers  corresponding 
to  the  data  shown  in  Figs.  (5.9)  and  (5.10).  It  may  be  supposed, 
therefore,  that  this  same  parameter,  t,  is  responsible  for  the  varia- 
tion in  the  results  of  other  workers. 

In  the  present  study,  no  independent  measurement  was  made  of  the 
actual  thickness  of  the  oxide  layer  on  any  of  the  samples.  Therefore, 
it  is  necessary  to  use  thickness  data  obtained  by  other  workers  for 
calibration.  Because  it  is  extremely  difficult  to  determine  the  actual 
thickness  of  the  oxide  layer  on  aluminum  substrates  and  because  the 
thickness  is  subject  to  variations  based  upon  the  sample  history,  there 
is  disagreement  in  the  available  literature.  Nevertheless,  certain 
limiting  values  can  be  presumed  with  reasonably  good  reliability  and 
such  limits  can  be  used,  in  conjunction  with  the  results  shown  in 
Fig.  (5.13),  to  deduce  a likely  thickness  scale  for  the  data  obtained 
in  this  study.  The  strategy  used  for  the  calibration  will  be  explained 
below. 

The  data  for  Sample  A and  Sample  B shown  in  Fig.  (5.7)  appear  to 
be  mutually  consistent  so  that  the  data  has  been  combined  into  a single 
set  for  the  succeeding  figures  and  for  conclusions,  such  as  Eqs.  (5.1) 
and  (5.2)  based  upon  those  figures.  The  fact  that  the  histories  of  the 
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two  samples  was  quite  different  will  now  be  used  to  determine  the 
thickness  scale.  Such  a procedure  can  only  be  justified  by  the  mutual 
consistency  evidenced  in  Fig.  (5.7):  a single  smooth  trajectory  can  be 
passed  through  the  data. 

Although  Sample  A and  Sample  B were  prepared  similarly.  Sample  A 
had  been  exposed  to  the  air  for  well  over  a year,  whereas  Sample  B was 
measured  within  four  hours  of  removal  of  the  collodion.  Thus  it  is 
reasonable  to  expect  that  the  layer  on  Sample  A had  reached  the  asymp- 
totic thickness  of  40  to  55A.  Although  this  asymptotic  thickness  can- 
not be  stated  with  great  precision,  many  workers  have  determined  it  to 
lie  in  this  range. [6  and  7]  The  layer  thickness  on  Sample  B,  on  the 
other  hand,  should  be  much  thinner  and  its  value  can  be  estimated  from 
the  results  obtained  by  Fane  and  Neal [6]  shown  in  Fig.  (5.1). 

The  results  shown  in  this  figure  allow  estimation  of  6A  as  a 
function  of  growth  time  of  the  oxide  layer,  and  show  a value  of 
6A  «100  minutes  of  angle  at  four  hours.  The  results  of  the  various 
studies  shown  in  Fig.  (5.12)  can  then  be  used  to  relate  6A  to  layer 
thickness.  Using  the  theoretical  thickness  scales  shown  on  the  figure, 
variations  in  A can  be  read  from  the  vertical  axis  corresponding  to 
thickness  values  from  50A  to  lOoA.  The  results  for  the  six  curves 
shown  are  (in  A/100  minutes  of  angle  in  A): 

10.62,  10.55,  10.75,  10.82,  11.11,  and  11.33. 

The  average  of  these  results  is 

fr-  10.9  A/100  minutes. 
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Thus  the  thickness  can  be  estimated  for  Sample  B: 

t £/  (|£)(5A)  d (10.9  S/100  min)(100  min)  10. 9S 
at  four  hours  of  growth  time. 

In  Fig.  (5.14)  the  theoretical  film  thickness  is  shown  versus  A 
for  the  two  extreme  index  values  marked  by  circles  on  Curve  B of 
Fig.  (5.13): 


ns  = 1.65(1  + j4. 70)  and 


ns  = 1.05(1  + j5.56) 


These  extreme  values  of  the  index  encompass  the  wide  variation  found  by 
previous  workers  in  the  field.  One  can  see,  directly  in  Fig.  (5.14), 
the  range  of  possible  value  of  t corresponding  to  a given  value  of  A. 
The  two  limiting  cases  discussed  above  for  Samples  A and  B before 
heating  are  also  marked  on  Fig,  (5.14)  as  rectangular  regions: 

Sample  A.  age:  over  a year 

thickness:  asymptotic,  40  to  55& 

A values:  five,  corresponding  to  Spots  1-5 

Sample  B.  age:  four  hours 

thickness:  less  than  12A,  approximately  10. 9^ 

A values:  two,  corresponding  to  Spots  a and  b. 

The  rectangular  regions  shown  in  Fig.  (5.14)  barely  intersect  the 
band  defined  by  the  index  extremes  shown.  The  third  curve  in  the 


figure,  corresponding  to 


n$  = 1.08(1  + j 5 . 555 ) 


(5.3) 


A 

Fig.  5.14  Theoretical  Film  Thic 
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was  selected  as  a compromise  between  the  index  data  found  in  this  study 
(and  by  previous  workers)  and  thickness  data  obtained  from  other 
workers.  The  point  from  Eq.  (5.3)  is  shown  as  a square  on  Curve  B in 
Fig.  (5.13). 

The  full  model  for  oxide  layers  on  pure  aluminum,  developed  in  the 
present  study,  is  summarized  by  Eqs.  (5.1),  (5.2),  and  (5.3).  The 
thickness  curve  corresponding  to  this  model  is  shown  on  the  \|>,  A plot 
of  Fig.  (5.15).  All  of  the  data  obtained  by  heating  Samples  A and  B 
is  shown  on  this  curve. 

The  justification  given  above  for  the  model  deduced  for  oxide 
layers  on  pure  aluminum  is  certainly  vulnerable  to  errors.  Based  upon 
the  intersection  of  the  rectangles  shown  in  Fig.  (5.14)  with  the  index 
band,  one  might  claim  an  uncertainty  in  t ranging  from 

t + 7.5X  to  t - 15A  . (5.4) 

On  the  other  hand,  one  might  wish  to  recognize  the  extremely  conserva- 
tive error  estimates  which  result  from  the  entire  band  of  index 
uncertainty.  Such  an  error  bound  would  vary  with  A and  would  range 
from 

t +_57'l\  at  A = 130° 

to  (5.5) 

O 

t at  A = 110° 


> 


Granted  the  possibility  of  absolute  error  in  the  construction  of 
the  model  of  this  study,  it  is  also  of  some  interest  to  apply  some 
significance  test  to  the  theory.  This  is,  after  all,  a rather  presump- 
tive theory,  presuming  that  a single  parameter  t,  defining  an  effective 
thickness  according  to  the  single  curve  shown  in  Fig.  (5.15)  can  use- 
fully characterize  the  \p,  A results  of  the  somewhat  rough  oxide  films 
grown  on  pure  aluminum  substrates. 

In  order  to  make  a crude  significance  test  of  the  model,  the 
following  procedure  was  followed,  leading  to  the  definition  of  a statis- 
tical parameter:  5tmocje-|-  Referring  to  Fig.  (5.15)  presume  that  some 
origin  is  selected  on  the  41,  A plot  at  iJj0  around  41°  or  42°  and  A0 
larger  than  the  value  of  the  oA  point  of  the  thickness  curve.  Call  the 
origin  selected  0.  The  exact  location  of  0 will  not  matter  to  this 
test  in  first  order. 

Now,  any  point  on  the  theory  curve  can  be  located  by  its  distance 

rc<tc> 

from  0,  where  tc  is  the  value  of  the  thickness  scale  reading  at  the 
selected  point  on  the  curve.  The  distance  rc  can  be  measured  in  any 
units,  say  inches,  because  its  only  relevant  property  is  that  rc(tc)  is 
a monotonically  increasing  function  of  tc.  Because  of  the  monotonicity, 
one  may  speak  instead  of  the  function  tc(rc).  Similarly,  the  ith  data 
point  on  Fig.  (5.15)  can  be  measured  at  a distance  r^  from  0.  A per- 
pendicular projection  of  the  data  point  onto  the  closest  point  of  the 
theory  curve  then  defines  a thickness  reading  ti  which  can  be 
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associated  with  ri  to  give  the  parametrization 

r^t.)  or  t^r.) 

for  each  data  point. 

The  curve  in  Fig.  (5.16)  shows  t (r  ) with  r in  inches,  measured 

L U V* 

from  an  origin,  0,  located  at  (ij;0,  A0)  = (41.5°,  155°).  Also  shown  in 
the  figure  are  all  of  the  points  t- (r.)  obtained  from  the  data  shown 
in  Fig.  (5.15).  The  data  points  in  Fig.  (5.16)  are  then  bounded  by 
drawing  two  curves  tu(rc)  and  T^(rc),  upper  and  lower  data  bounds, 
respectively.  In  order  now  to  measure  the  significance  of  the  theory 
with  respect  to  the  data,  the  parameter  rc  can  be  eliminated  graphical- 
ly and  one  can  obtain  ty  and  t^  as  functions  of  tc-  Then  the 
difference 

6ttheory^V  E W " W 

may  be  obtained  as  a function  of  tc.  This  difference  is  shown  in 

O 

Fig.  (5.17),  expressed  in  A.  It  must  be  stressed  that  the  quantity 
defined  in  Eq.  (5.6)  is  not  an  error  in  the  usual  sense  of  error 
analysis,  but  is,  rather,  a measure  of  the  relevance  of  the  theory  to 
the  data  supporting  it.  This  theory  significance  measure  is  expressed 

O 

in  A simply  to  offer  some  perspective  to  the  measure.  The  results 
shown  in  Fig.  (5.17)  have  only  the  qualitative  interpretation  that  the 
theory  is  reasonably  well  matched  to  the  data.  Another  theory  curve 
tested  in  the  same  way  against  the  same  data  might  have  revealed 
considerably  larger  values  of  6 ttheory  and,  therefore,  might  be  judged 


to  be  a less  significant  theory.  Also,  the  width  of  the  “theory  error 
band"  shown  in  Fig.  (5.17)  might  indicate  a systematic  increase  as, 
say,  tc  increased.  Such  behavior  would  offer  the  qualitative  sugges- 
tion that  the  theory  begins  to  fail  for  large  tc- 

As  one  examines  Fig.  (5.15),  it  may  be  noticed  that  there  is  a 
systematic  error  in  the  fit  to  the  data:  for  large  A the  data  tends  to 
be  to  the  right  of  the  curve;  for  small  A the  data  tends  to  be  to  the 
left  of  the  curve.  It  is  obvious,  therefore,  that  a better  curve  fit 
could  be  achieved.  The  St^^  test  proposed  above  would,  however, 
give  about  the  same  results  for  a better  curve  fit.  Goodness  of  data 
fit  is  measured  by  quantities  such  as  RMS  deviation  of  the  data  from 
the  curve.  The  ^t^eo  test  is  intended  to  try  to  measure  another 
feature  of  the  theory:  is  it  a useful  way  to  parametrize  the  data? 
Although  a better  curve  fit  is  obviously  possible,  there  is  not  a large 
enough  data  base  in  the  present  study  to  justify  additional  tuning.  It 
does  appear,  however,  that  the  strategy  of  the  theory  is  meaningful. 

Because  of  the  reasonably  complex  method  used  to  determine  the 
various  effective  parameters  which  characterize  the  media,  it  is  useful 
to  note  the  sensitivity  of  the  thickness  parameter  to  the  various 
parameters.  In  Figs.  (5.18),  (5.19),  and  (5.20)  the  model  curve, 
corresponding  to  Eqs.  (5.1)  to  (5.3),  is  shown  together  with  nearby 
curves  that  would  result  if  one  of  the  media  parameters  were  to  be 
charged  slightly.  In  Fig.  (5.18)  variations  in  Kf  are  illustrated  by 
modifying  Eq.  (5.1)  to: 

Kf  = (0.01)^  ± 0.01 
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to  obtain  the  two  additional  curves. 

In  Fig.  (5.19)  the  nominal  value  of  5.555  for  K$  is  varied  to  6.0 
and  5.0  for  comparison.  Fig.  (5.20)  shows  four  additional  values  for 
n$  along  with  the  nominal  value  of  1.08.  The  parameter  n^  is  not 
illustrated  here,  because  it  was  arbitrarily  locked  at  1.65  for  the 
present  study.  Although  A is  slightly  sensitive  to  n^,  the  parameter 
n^  enters  in  such  a way  that  the  conclusions  of  this  study  are 
relatively  insensitive  to  this  parameter. 

The  sensitivity  results  illustrated  in  Figs.  (5.18)  to  (5.20)  may 
be  suimiarized  in  terms  of  the  partial  derivatives  of  A,  ij/,  and  t listed 
in  Tables  (5.2),  (5.3),  and  (5.4).  Variation  of  the  incidence  angle  0^, 
about  the  nominal  value  of  70°,  is  illustrated  in  Fig.  (5.21)  by 
nominal  media  parameter  curves  of 

6.  * 70°  ± 0.2°  . 

Partial  derivatives  of  A and  ip  with  respect  to  0.,  evaluated  from 
these  results,  are  listed  in  Table  (5.5).  Corresponding  partial 
derivatives  of  t with  respect  to  0^  may  be  obtained  by  combing  resuHs 
from  Table  (5.3)  and  Table  (5.6): 

6t  _ frt  , frt  6ip 
60.  Sip  60 i 6A  66i 

To  use  the  model  shown  in  Fig.  (5.15)  a data  point  is  projected  in 
a perpendicular  fashion  to  the  closest  point  on  the  theory  curve  and 
the  corresponding  value  of  t is  then  obtained  from  the  scale  on  the 
theory  curve.  In  a sense,  errors  due  to  media  parameters,  alignment 
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TABLE  5.3 

SENSITIVITY  OF  THICKNESS  TO  A AND  ip 


9 

Range  of  t (A) 

fit  , JL 
fi^  ^ 

fit  {A 
6A  M° 

0 - 

50 

3.02 

6.35 

« 

50  - 

100 

1.43 

6.78 

100  - 

150 

0 

7.27 

150  - 

200 

1.96 

7.84 

9 

200  - 

250 

5.55 

8.89 

250  - 

300 

6.83 

9.76 

300  - 

350 

9.72 

11.11 

9 


9 


9 


9 
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TABLE  5.4 

SENSITIVITY  OF  A TO  K$,  n$  AND  Kf 


Value„of  t 6A  6A  6A  6t 

(A)  6Ks  6n$  6Kf  6 A 


0 

6.84 

37.32 

0 

6.35 

50 

6.48 

35.34 

.4 

6.57 

100 

6.10 

33.18 

1.2 

7.03 

150 

5.69 

30.94 

2.3 

7.56 

200 

5.28 

28.66 

3.75 

8.37 

250 

4.87 

26.38 

5.3 

9.33 

300 

4.46 

7.05 

10.44 

350 

4.06 

9.25 
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TABLE  5.6 

SENSITIVITY  OF  THICKNESS  TO  6i 


Value  of  t(A) 

6iJ> 

6A 

60. 

60i 

0 

.15 

2.47 

50 

.12 

2.78 

100 

.10 

3.03 

150 

.078 

3.21 

200 

.055 

3.34 

250 

.0375 

3.42 

300 

.0175 

3.47 

350 

.0025 

3.49 

errors,  and  the  random  errors  arising  in  the  taking  of  the  data  have 
already  been  incorporated  into  the  model  developed  in  Eqs.  (5.1)  to 
(5.3).  Thus,  it  is  not  possible  to  untangle  expected  errors  more 
accurately  than  the  large  and  small  error  estimates  made  in  Eqs.  (5.5) 
and  (5.4),  respectively. 


SECTION  VI 


A STUDY  OF  OXIDE  FILMS  ON  ALCLAD  ALLOY 

Once  the  procedure  for  the  ideal,  pure  aluminum  surfaces, 
described  in  Section  V,  was  developed,  measurements  were  attempted  for 
the  more  difficult  case  of  alclad  aluminum  alloy.  The  intention  was  to 
attempt  to  develop  a method  for  using  ellipsometric  measurements  to 
characterize  such  rough  surfaces.  It  is  of  considerable  interest  to 
determine  the  impact  of  surface  roughness  on  the  values  of  n$  and  n^  as 
determined  by  ellipsometry  and  to  attempt  to  determine  the  effective 
thickness  of  the  oxide  layer. 

Alclad  2024,  the  best  known  and  most  widely  used  aircraft  alloy, 
was  used  for  the  ellipsometric  measurements.  This  alclad  has  a core  of 
alloy  2024  with  copper  as  the  principal  alloying  element.  The  elements 
in  this  alloy  are: 

.50%  Silicon 
.5%  Iron 

3.8  - 4.9%  Copper 
.30  - .9%  Manganese 
.1%  Chromium 
.25%  Zinc 
.15%  Titanium 
.15%  Other 
remainder  Aluminum 

The  core  is  covered  on  both  sides  with  a cladding  of  alloy  1230.  The 
minimum  thickness  of  the  cladding  is  2%  of  the  total  thickness  per  side. 
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The  cladding  is  composed  of  a minimum  of  99.3%  aluminum  with  other 
elements  being 
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.7%  Silicon  and  Iron 
.10%  Copper 
.05%  Manganese 
.1%  Zinc 

.05%  Other  elements 


The  thickness  of  the  alclad  used  was  approximately  .04  inches.  The 
alclad  is  solution  heat-treated  at  493°C  and  then  worked  cold  to 
improve  strength. 

As  had  been  done  with  the  pure  aluminum  samples,  it  was  decided  to 
heat  the  alclad  samples  over  and  over,  determining  the  parameters  and 
A after  each  heating.  In  this  way  a data  trajectory  can  be  established 
on  a i|i,  A plot  and  selection  of  effective  parameters  to  fit  the  data 
can  be  attempted.  The  melting  range  of  alclad  is  between  502°  and 
638°C  so  that  a heating  temperature  of  350°C  was  again  selected.  After 
initial  warmup,  the  samples  were  held  at  350°C  for  two  hours  and  then 
allowed  to  cool. 

On  the  first  sample.  Sample  1,  four  spots,  a,  b,  c,  and  d,  were 
marked.  The  measurements  for  \p  and  A were  made  at  these  points,  then 
the  sample  was  heated  four  times  with  measurements  made  after  each 
heating.  Data  points  for  spots  a,  b,  c,  and  d are  given  subscripts 
0,  1,  2,  3,  4,  denoting  the  number  of  times  the  sample  had  been  heated 
before  the  measurement.  After  the  first  heating,  the  measured  value 
for  A increased  for  three  of  the  points.  Further  heating  did  decrease 
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I the  measured  A as  expected  from  experience  with  the  optically  flat 

aluminum  samples.  As  shown  in  Fig.  (6.1)  the  data  from  the  measure- 
ments are  scattered  over  a much  wider  range  than  was  found  for  aluminum 

> samples  in  the  previous  section,  although,  as  the  samples  are  heated 
more  times,  they  tend  to  develop  parallel  trajectories  of  successively 
decreasing  A. 

♦ Careful  examination  of  the  alclad  shows  very  shallow  grooves  in 
the  surface:  rolling  grain  marks  left  from  the  rollers  used  in  manu- 
facturing. Measurements  made  with  these  lines  vertical,  designated  as 

♦ a 90°  sample  angle,  result  in  the  light  shining  across  these  grooves 
and  the  surface,  perhaps,  appears  to  the  light  to  be  rougher.  Measure- 
ments made  with  these  grooves  horizontal,  designated  as  0°  sample  angle, 

> effectively  appear  to  be  a smoother  surface.  In  a sense,  the  grooves 
generate  a family  of  parallel  smoother  surfaces.  Measurements  were 
taken  from  two  different  samples  to  find  if,  indeed,  the  sample  angle 

♦ has  any  effect  on  the  values  of  and  A.  The  measurements  from 
Sample  2,  Fig.  (6.2)  were  taken  at  two  spots,  a and  b.  Measurements  at 
each  spot  were  taken  at  both  0°  and  90°  sample  angles  before  heating 

♦ and  after  three  successive  heatings.  Again,  subscripts  are  used  to 
indicate  the  number  of  heating  cycles.  The  data  for  the  smoother  sur- 
face, 0°  sample  angle,  has  larger  values  of  i/>,  although  the  trajectory 

♦ for  a given  point  on  the  sample  follows  a similar  path  for  both  the 
0°  and  90°  angles. 

» 
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Measurements  from  a third  piece  of  alclad.  Sample  3,  gave  results 
similar  to  those  for  Sample  2.  These  data  are  shown  in  Fig.  (6.3).  A 
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Fig.  6.2  Data  from  Alclad  Sample  2 Taken  at  Sample  Angles  of 
Both  0°  and  90°  with  Trajectories  of  Increasing 
Thickness 


* 

Fig.  6.3  Data  from  Alclad  Sample  3 Taken  at  Sample  Angles  of 
both  0°  and  90°  with  Trajectories  of  Increasing 
Thickness 
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shift  of  both  the  90°  sample  angle  trajectory  and  the  0°  sample  angle 
trajectory  to  the  right  (larger  values  of  \J»)  for  Sample  3 with  respect 
to  Sample  2 can  be  observed  by  comparing  Figs.  (6.2)  and  (6.3).  From 
the  work  of  Section  V,  the  shift  to  larger  \p  is  known  to  be  correlated 
with  a smaller  effective  roughness.  As  a matter  of  fact,  the  surface 
of  Sample  3 looked  smoother  and  shinier  than  the  surface  of  Sample  2. 

At  a single  point  on  a sample  the  thickness  is  obviously  the  same 
whether  measured  at  a sample  angle  of  0°  or  at  a sample  angle  of  90°. 
Using  the  measurements  taken  at  angles  of  0 and  90°  for  each  point  on 
the  two  samples,  lines  of  constant  thickness  can  be  drawn.  The  results 
are  shown  in  Figs.  (6.4)  and  (6.5).  Using  the  trajectories  formed  from 
the  data  taken  after  consecutive  heatings,  shown  in  Figs.  (6.2)  and 
(6.3),  and  the  lines  of  constant  thickness,  from  Figs.  (6.4)  and  (6.5), 
an  attempt  was  made  to  fit  a theory  curve  to  the  measured  data. 

The  first  attempts  in  fitting  a curve  to  the  alclad  data  involved 
the  variation  of  only  one  parameter  from  the  set  used  for  the  optically 
flat  samples  of  aluminum.  The  first  parameter  varied  was  n$  while 
K$  = 5.555,  n^  = 1.65  and  K ^ = (.01)^  were  kept  fixed.  The  value  of  ns 
was  varied  from  0.6  to  1.08  and  the  curves  of  constant  thickness  for 
this  variation  of  ns  are  shown  in  Fig.  (6.6).  Also  shown  in  Fig.  (6.6) 
are  data  measured  at  spot  a of  alclad  Sample  3.  Comparison  of  the  cal- 
culated curves  with  the  measured  data  from  Sample  3 shows  that  the 
slope  of  the  calculated  constant  thickness  lines  are  much  steeper  than 
the  constant  thickness  lines  for  the  measured  data. 

Next,  the  parameter  K was  varied  while  n$  was  returned  to  1.08. 
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The  range  for  which  K$  varied  was  from  1.45  to  1.60.  The  curves 
produced  from  this  variation  are  shown  in  Fig.  (6.7)  along  with  the 
data  from  spot  a of  Sample  3.  Comparison  with  the  data  of  the  constant 
thickness  lines  produced  from  variations  of  K%  shows  that  these  theore- 
tical curves  are  also  too  steep  to  fit  the  data. 

A combination  using  both  K$  and  n$  was  then  tried.  By  judging  the 
behavior  of  the  curves  in  Figs.  (6.6)  and  (6.7)  and  by  experimenting 
with  the  parameter  dependence  of  the  curves,  the  following  model  was 
devised: 

K = 9.75  - 4. 77n  + 0.505n*  . (6.1) 

Curves  based  upon  this  relationship  are  shown  in  Fig.  (6.8).  Again 
comparison  can  be  made  to  the  data  taken  from  spot  a of  Sample  3.  The 
slope  of  the  constant  thickness  lines  for  the  theory  curves  and  the 
data  do  agree,  but  the  trajectories  of  increasing  oxide  thickness  do 
not  follow  the  lines  of  constant  n$. 

In  Section  V,  a change  in  K ^ was  used  to  correct  for  variation  be- 
tween the  theory  and  the  data.  The  change  was  made  to  correct  for 
surface  roughness  and  was  varied  as  a function  of  thickness: 

Kf=(°.01)^)  • (6.2) 

This  same  form  of  Kf  as  a function  of  t was  used  for  the  curves  in 
Fig.  (6.8). 

It  might  be  assumed  that  the  roughness  of  the  oxide  layer  growth 
responds  to  increases  in  the  substrate  roughness.  A way,  therefore. 
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to  improve  the  curve  fit  might  be  to  increase  K ^ in  response,  not  only 
to  t,  but,  also,  in  response  to  substrate  roughness  as  measured  by  K$. 
However,  since  Eq.  (6.1)  models  K$  as  a function  of  n$,  it  is  equiva- 
lent to  assign  a dependence  on  n$  as  well  as  t.  In  this  spirit,  the 
following  model  was  tried: 

Kf=(°.01)^ns  • (6.3) 

The  curves  displayed  in  Fig.  (6.9)  correspond  to  the  same  case 
shown  in  Fig.  (6.8)  except  that  Eq.  (6.3)  is  used  to  describe  Kf.  The 
data  from  both  spots  on  Sample  3 are  shown  for  comparison.  The  data 
measured  from  Sample  2 are  shown  with  the  same  theory  curves  in 
Fig.  (6.10).  The  curves  shown  in  Figs.  (6.9)  and  (6.10)  are  obtained 
by  letting  ns  vary  from  1.08  to  2.2  while  K$  and  vary  as  functions 
of  n$  according  to  Eqs.  (6.1)  and  (6.3).  Fig.  (6.11)  shows  the  func- 
tions and  defined  by  Eqs.  (6.1)  and  (6.3). 

Again  reference  is  made  to  the  work  of  Fenstermaker  and 
McCracken [9]  showing  the  theoretical  change  of  n$  and  K$  as  a surface 
of  gold  or  silver  increases  in  roughness.  Their  results,  reproduced  in 
Fig.  (6.12)  shows  an  increase  in  n$  and  a decrease  in  Ks  as  the  surface 
roughness  increased.  As  evidenced  in  Fig.  (6.11),  this  is  the  type  and 
direction  of  the  variations  used  in  fitting  the  present  data  of  alclad. 

Now  that  a model  for  the  parameter  behavior  for  oxide  films  on 
alclad  has  been  selected  to  fit  the  data  trajectories,  certain  obser- 
vations can  be  made.  First,  it  does  appear  that  the  thickness  of  oxide 
layers  on  aluminum  can  be  measured  with  ell ipsometry,  provided  that 
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Fig.  6.10  Lines  of  Constant  Thickness  Derived  from  Variations 
in  n$  where  Ks  and  Kf  are  Functions  of  ns  Compared 
with  Data  from  Sample  2 


I 


124 


» 


» 


♦ 


> 


I 


» 


♦ 


I 


» 


» 


Surface  Roughness,  A 


Fig.  6.12  Apparent  Values  of  n and  K for  Rough  Surfaces 
Derived  using  Triangular,  Square  and  Pyramid 
Models 
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attention  is  given  to  the  sample  angle.  As  is  evident  in  Fig.  (6.9), 
for  Sample  2 the  effective  n$  for  a 90°  sample  angle  differs  from  the 
effective  ns  for  a 0°  sample  angle  by  about  0.4.  A similar  difference 
is  found  for  Sample  3 in  Fig.  (6.10).  One  notices  a larger  average 
value  of  n$  for  Sample  2,  the  one  appearing  rougher  to  the  eye. 

As  originally  anticipated,  the  sample  angle  of  90°  is  shown,  by 
Figs.  (6.9)  and  (6.10),  to  correspond  to  a greater  roughness  assign- 
ment (larger  value  of  n$)  than  the  sample  angle  of  0°.  This  conclusion 
was  tested  at  intermediate  sample  angles  for  both  Samples  2 and  3.  The 
results,  shown  in  Fig.  (6.13),  indicate  that  one  can  almost  use  the 
ellipsometrically  determined  value  of  n$  to  detect  the  sample  orienta- 
tion angle,  if  the  average  value  of  n$  is  known  for  the  sample. 

The  results  in  Fig.  (6.13)  were  taken  after  Samples  2 and  3 had 
been  heated  thrre  times  for  a total  of  six  hours.  The  same  effect  was 
examined  on  a fresh,  unheated,  sample  of  alclad  and  the  results  are 
shown  for  two  points  on  the  sample  in  Fig.  (6.14). 

In  conclusion,  it  appears  that  the  method  employed  in  this  study, 
using  data  trajectories  obtained  by  repeatedly  heating  samples  to 
deduce  parameter  interrelationships,  shows  great  promise  for  ellipso- 
metric  work  for  rough  films  on  both  smooth  and  rough  substrates.  The 
measurements  do  indicate  that  index  parameters  change  in  an  orderly 
way  that  can  be  predicted  as  functions  o^  thickness  and  surface  rough- 
ness. Additional  work  will  be  required  to  test  the  specific  parameter 
relationships  found  in  this  work. 
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